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ABSTRACT 

We present a spectroscopic survey of 318 faint iJZ ~ 27, L ~ O.IL*), Lya-emission-selected galaxies 
(LAEs) in regions centered on the positions of hyperluminous QSOs (HLQSOs) at 2.5 < z < 3. A 
sample of 32 LAEs with rest-frame optical emission line spectra from Aecfc/MOSFIRE are used to 
interpret the LAE spectra in the context of their systemic redshifts. The fields are part of the Keck 
Baryonic Structure Survey (KBSS), which includes substantial ancillary multi-wavelength imaging 
from both the ground and space. From a quantitative analysis of the diverse Lya spectral morpholo¬ 
gies, including line widths, asymmetries, and multi-peaked profiles, we find that peak widths and 
separations are typically smaller than among samples of more luminous continuum-selected galaxies 
(LBGs) at similar redshifts. We find tentative evidence for an association between Lya spectral mor¬ 
phology and external illumination by the nearby HLQSO. Using the MOSFIRE subsample, we find 
that the peak of the resolved (i? « 1300) Lya line is shifted by -1-200 km s“^ with respect to systemic 
across a diverse set of galaxies including both LAEs and LBGs. We also find a small number of 
objects with significantly blueshifted Lya emission, a potential indicator of accreting gas. The Lya- 
to-Ha line ratios measured for the MOSFIRE subset suggest that the LAEs in this sample have Lya 
escape fractions /esc,Lya ~ 30%, significantly higher than typical LBG samples. Using redshifts cali¬ 
brated by our MOSFIRE sample, we construct composite LAE spectra, finding the first evidence for 
metal-enriched outflows in such intrinsically-faint high-redshift galaxies. These outflows have smaller 
continuum covering fractions (/c « 0.3) and velocities (uave ~ 100 — 200 km s“^, Umax ~ 500 km s“^) 
than those associated with typical LBGs, suggesting that gas covering fraction is a likely driver of the 
high Lya and Ly-continuum escape fractions of LAEs with respect to LBGs. Our results suggest a 
similar scaling of outflow velocity with star formation rate as is observed at lower redshifts (uoutflow 
SFR°-^®) and indicate that a substantial fraction of gas is ejected with v > Vesc- Further observations, 
including deep spectroscopy in the observed near-IR, will further probe the evolution and enrichment 
of these galaxies in the context of their gaseous environments. 

Subject headings: galaxies: formation — galaxies: high-redshift — galaxies: intergalactic medium 


1. INTRODUCTION 


Stellar feedback plays a crucial role in galaxy forma¬ 
tion. Massive stars enrich, heat, ionize, and expel gas 
throughout their lives and deaths, and understanding 
these processes is key to modeling the observational char¬ 
acteristics of galaxies at all redshifts. These character¬ 
istics include t he low instantan eous efficiency of star- 


and baryon content of galaxies (e.g., 

White & Frenk 1991 

Keres et al. 2009 

Hopkins et al. 

2U14P, and the pres- 

ence and evolution ot the mass-metaliicitv relationship 

of galaxies (e.g., 

Erb et al. 

2006a 

Andrews & Martini 
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2013| Steidel et al. 20141. Each of these observations 
requires that the majority of gas associated with galaxy 
halos is prohibited from forming stars, particularly in the 
lowest-mass galaxies. There are many potential causes 
for this inefficiency: the temperature, ionization, and/or 
kinematics of the gas may suppress its conversion into 
stars. Simulations can reproduce many characteristics of 
the stellar populations of galaxies simply by artificially 
delaying gas cooling or implementing other ad hoc solu¬ 
tions specific to these discrepancies; these naturally lack 
the predictive power of a feedback model based on the 
phys ics of star formation and supernovae thems elves (see 


e.g., |Hummels et al.||2013t |Hopkins et al. 
In particular, mahy simulations strugg. 


2014). 


e to reproduce 
properties of the circumgalactic medium (GGM; e.g.. 


Rudie et al. 2012) of galaxies, which must be closely 
linked to the same feedback processes that shape their 
stellar populations. Studies of the gas in and around 
galaxy halos conclusively show that galaxies enrich their 
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1 Aguirre et al.|2001 

Pettini et al.|2003 

Adelberger et al. 


2UU5t Steidel e 

al.|20101|Martin et al.|2010 

'Lurner 

et al. 

2014 

, and it is expected that the chemistry and 


enriched gaseous outflows seen ubiqit ously in the spectra 


of high-redshlft star-forming galaxies dSteidel et al. 

1996 

Pettini et al. 

20011 Shapley et al. 120031 

Jones et al. 


Martin et al. 

2012). 


after referred to as LBGs), these outflows produce spec¬ 
tral signatures In both absorption and emission, allowing 
the characterization of the chemistry and kinem atics of 
the outflows for gravitatlonally-lensed s ources (j Pettim 


et al. I 2 OOOI I 2 OO 2 ), or sta cks of galaxies (Shapley et al. 
2U03 [steidel et al'||2010| . Lyman-a emitters (LAEs) as 


a population a re found to be youn ger, fainter, and less 
reddened (e.g., Gawlser et al. [20061 than LBGs on aver- 


ageHso LAE surveys provide a complementary sample to 
LBGs for characterizing the effects of these gas flows (and 
stellar feedback in general) on galaxies over a broad range 
of ages and masses. Gompared to most LBGs, LAEs 
may also “leak” a relatively large fraction of the Lyman- 
contlnuum (LyG) ph otons produced In their embedded 
star-f orming regions (Nestor et al.||2013 Mostardi et al. 
20131, which may be linked to the observed association 
between higher Lya equivalent widt hs and lower neutral 
gas covering fract ions among LBGs ( Shapley et al. [20(33 
Jones et al.|2013|). Therefore, LAEs may be most analo¬ 
gous to those gala xies that dominate the relonlzation of 
the universe (e.g., Robertson et al. 2013 Dressier et al. 
20T4|. 


Most Lya-selected galaxies are too faint for contin¬ 
uum absorption features to be easily observed. At 


low redshifts, the Lya Referen ce Sample (LARS; Hayes 


et al. 


2013 


Ostlin et al. 


20141 includes 14 local galax¬ 


ies with HS'l' imaging and spectroscop y of their Lya 


emission. 


Rivera-Thorsen et al. (20151 analyze metal 


absorption and Lya emission troin these galaxies, find¬ 
ing that gaseous outflows and the covering fractions 
of absorbers play a strong role i n modulating Lya es - 


Wofford et al. 


2013). 


cape (see also similar results by 
The few measurements of interstellar absorption in high- 
reds hift LAE spectra come f r om especially bright exam¬ 
ples (Hashlmoto et al. 2013| |Shlbuya et al. 2014), with 
continuum magnitudes (i^ ~ 24, L ~ L*j bright enough 
to fall Into typical samples of LBGsj^ However, it is the 
faint LAEs that provide the most interesting environ¬ 
ment for studying the kinematics of gas and galaxies. As¬ 
suming that the feedback processes ubiquitous in LBGs 
extend to the earliest and faintest galaxies, then the effect 
of even faint galaxies on their surrounding gaseous envi¬ 
ronments must be taken into account when modeling the 
enrichment of the circ umgalactic and inter galactic media 
(GGM and IGM; e.g., Steidel et ah 2010) and its trans¬ 


parency to ionizing photons (e.g., Schenker et al. 2013 
Rudie et al. 2013). Furthermore, the low masses anc 


^ Note, however, that all LBGs emit at least low-level Lya emis¬ 
sion, and all LAEs would display Ly man breaks in suffic iently sen¬ 
sitive broadband images or spectra llSteidel et al.|2011||. _ 

® For example, 25% of the 811 “LBGs" m |bhapley et al.|j2003| 
have WLyo, > 20A and spectroscopic properties similar to these 
bright-LAE samples. 


star formation rates of galaxies selected as LAEs may 
enable t he detection of the gas Infall predicted by s lmu- 


lations (Birnboim & Dekel 2003 Keres et al. 2005 and 
others) whereas in LBGs with stronger star-Tormation 
rates, gaseous outflows may be more likely to dominate 
over infall In determining the state of circumgalactic gas. 

Regions of the sky around bright QSOs provide par¬ 
ticularly advantageous environments for LAE surveys. 
Hyp erlumlnous QSOs inhabi t known galaxy overdensi¬ 


ties (Trainor & Steidel 2012) and thereby yield a large 
number of LAEs within a single telescope pointing and 
narrowband-selected redshlft Interval. Additionally, the 
Ionizing fields of the QSOs can boost the Lya emission 
of faint, gas-rich galaxies through Lya “fluorescence” 


1 Hogan & Weymann 

1987 Adelberger et al. 

2006 Can- 

talupo et al. 120071120 

[2 IKollmeier et al.|2010 

I'Lrainor &: 


ethcient tor the selection and follow-up of faint LAEs. 
Fluorescent Lya emission can allow the chara cterization 


of the QSO radiative geometry and history (Trainor & 


Steidel 2013), but it also provides a means of selecting 


galaxies irrespective of their intrinsic stellar luminosities. 
LAE s selected In these regio ns may include “dark galax¬ 
ies” (Gantalupo et al. 2012), halos of gas and dark mat¬ 


ter in which only negligible star formation has occurred. 
Whether or not such “dark” objects are observable, the 
fluorescent boost of nearby QSOs may allow us to probe 
an even greater range of galaxy properties than typical 
LAE surveys alone. 

In this paper, we present rest-frame UV spectroscopy 
from a large (68,000 cMpc^) narrowband survey con¬ 
ducted In eight QSO fields with 2.55 < zqso < 2.85. 
All eight fields are part of the Ke ck Baryonic Structure 


Survey (KBSS; Rudie et al. 2012), and this Lya survey 
(designated KBSS-Lya) is complementary to the original 
KBSS for probing the interactions of QSOs, galaxies, and 
gas at the peak redshifts of star formation. The authors 
note that the effect of the bright QSOs in the KBSS-Lya 
fields on their nearby galaxies Is not generally known, but 
the QSO-lnduced contribution to the Lya emission is ex¬ 
pected to be independent of the intrinsic star-formation 
properties of these galaxies. In this paper, we generally 
assume that the kinematic and transmissive properties of 
the KBSS-Lya LAEs are most directly related to their 
internal star formation and stellar populations. However, 
the results of our analysis that may depend strongly on 
QSO Illumination and proximity are noted explicitly. Fu¬ 
ture work will consider the effect of the QSQ on associ¬ 
ated galaxies in greater detail. 

The spectroscopic component of the KBSS-Lya sur¬ 
vey Is described in Sec. In Sec. we present the 
properties of the Lya emission lines, iiiGuding analysis of 
the Lya escape fraction, spectral morphology, and kine¬ 
matics in relation to comparison samples of continuum- 
bright galaxies. In Sec. we present measurements of 
metal absorption lines In the stacked continuum spec¬ 
tra of the rest-UV spectroscopic sample and characterize 
the ionization, covering fraction, and kinematics of star- 
formation driven outflows in the LAEs. A brief discus¬ 
sion of these results in the context of the physical prop¬ 
erties of the galaxies and previous work is given in Sec.[^ 
and our results are summarized in Sec. Physical quan¬ 
tities are given in observable units when possible, but a 
AGDM universe with (fim, Hq) = (0.3, 0.7, 70 km 
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TABLE 1 

Lyo-Emitter Field Descriptions 


QSO Field 

(±0.001) 

NB 

Filter 

BB 

Filter 

^spec 

QOlOO+13 (PHL957) 

2.721 

NB4535 

B+G 

20 

HS0105-I-1619 

2.652 

NB4430 

B 

22 

Q0142-10 (UM673a) 

2.743 

NB4535 

B+G 

22 

Q1009+29 (CSO 38) 

2.652 

NB4430 

B 

35 

HS1442-I-2931 

2.660 

NB4430 

B 

41 

HS1549-I-1919 

2.843 

NB4670 

G 

95 

HS1700-I-6416 

2.751 

NB4535 

B+G 

20 

Q2343+12 

2.573 

NB4325 

B 

63 


s is assumed when necessary. 


2. OBSERVATIONS 
2.1. LAE Sample 


The full parent sample of ^800 photometrically- 
identified LAEs will be described separately (R. Trainor 
et ah, in p rep.); this LAE sample w as also used and dis¬ 
cussed by |Trainor & Steidel| (|2013|) . LAEs were identi¬ 
fied on the basis of paired narrowband and broadband 
observations (from Keck 1/LRIS-B) using a custom set 
of narrowband filters. A total of four narrowband filters 
were used, each with a bandpass constructed to select the 
wavelength of Lya at one or more QSO redshifts (Eig.[^. 
Objects were identified and photometry was performed 
using SExtractor in one- and two-image modes to ex¬ 
tract narrowband and broadband magnitudes (uinb and 
tobb) sampling Lya. Hereafter, tobb refers to the mea¬ 
sured LAE magnitude in the B and/or G hlter used to 
estimate the continuum emission near Lya, whichever is 
better matched at a given redshift (Table . Because 
the B/G filter bandpasses include the Lya emission line, 
we use 77.-band measurements to sample the continuum 
without Lya contamination. We conducted follow-up 
spectroscopy of the subset of objects identified to have a 
narrowband color excess ttibb — wnb >0.6, which cor¬ 
responds to a limit in rest-frame Lya equivalent width 


WAya ^ 20A. The success of our spectroscopic follow-up 
observations dropped significantly for tonb > 26.5, so 
this limit (approximately the 3-(t depth of our narrow- 
band observations) was adopted as the faint cut-off for 
the photometric and spectroscopic samples. This mag¬ 
nitude limit corresponds to an integrated narrowband 
flux ELya « 1 X 10“^^ erg s“^ cm“^ for a continuum- 
free Lya emission line, or ELyc « 6 x 10“^® erg s“^ 
cm“^ for an emission line with our median value of 
(bbLya) ~ 40A. The apparent magnitude distribution of 
both the parent sample and the spectroscopic subsample 
of LAEs discussed in this paper are given in Eig.[^ Both 
the spectroscopic and photometric samples have median 


(”^ab(6930A)) « 27. A ccording to the gala xy luminosity 
functions measured by Reddy et al. (2008), this contin¬ 
uum magnitude corresponds to L ~ 0.1L/ at z ~ 2.7. 
Details of the QSO fields and the spectroscopic sample 
are given in Table 


2.2. UV Spectroscopic Observations 

Rest-UV spectra were obtained with Keck 1 /LRIS-B in 
the multislit mode using the 600/4000 grism and 560nm 
dichroic. All objects were observed with 1"2 wide slits, 
but the spectral resolution is limited by the 0'.'6—0"8 see¬ 


ing disk for these (typically spatially unresolved) LAEs. 
The effective resolution near the observed Lya wave¬ 
length for these spectra is i? ~ 1300, corresponding to 
a velocity resolution (T„ ^ 100 km s“^. Spectroscopic 
observations were performed in sets of 3—4 1800 s expo¬ 
sures, for a total of 1.5—2 hours of total integration time 
per object. By the completion of the survey, we had ob¬ 
tained spectra for 415 of the 749 photometric LAE can¬ 
didates, with some candidates observed on as many as 
three different masks. The observations were conducted 
and reduced in essentially the same manner as those of 
the KBSS and precedin g surveys, a detailed des cription 
of which can be found in|Steidel et al.l (|2003|) andlSteidell 


et al.| (20041. Specifically, masks were construc t ed an d 


targets assigned as described by |Steidel et al. (|2003|) 
but with a minimum slit length of 11" in order to pro¬ 
vide better background subtraction for continuum-faint 
{TZ ~ 27) sources. The data were reduced (including flat- 
fielding, cosmic-ray rejection, background subtraction, 
flexure compensation, and wavelength/flux calibration) 


using a cust om set of IRAE scripts as described in Steidel 


et al. (2003). The two-dimensional spectra were rectihed 


using hat-rield observations, and one-dimensional spec¬ 
tra were then extracted using the IRAF task apall with 
10 pixel apertures (1"35, ~2x the seeing FWHM). We 
detect little or no continuum emission in individual spec¬ 
tra, so all pixels in the aperture received equal weight¬ 
ing during the extraction and the object trace was not 
allowed to vary with wavelength along the rectified 2D 
spectrum. Note that the Keck 1 Cassegrain Atmospheric 
Dispersion Corrector (ADC) removes the wavelength- 
dependent shifts in object position due to differential at¬ 
mospheric refraction. The two-dimensional spectra were 
binned in two-pixel blocks along the wavelength axis dur¬ 
ing readout, and the one-dimensional extracted spectra 
were resampled to the median resulting wavelength scale 


of the observations, 1.26 A/pix. 

The resulting one-dimensional spectra were then sub¬ 
jected to an automated line-detection algorithm de¬ 
scribed here. First, each spectrum was smoothed with 
a kernel corresponding to the instrumental resolution 
(^3.5A FWHM). A detection region of the spectrum was 
then isolated: in order to ensure that spectroscopically- 
detected objects corresponded to their NB-selected coun¬ 
terparts, the algorithm only located emission lines that 
would fall between the 10% power points of the corre¬ 
sponding NB filter used to select them. The highest peak 
in this detection region was then used to define the first 
guess at the emission line wavelength. The significance 
of the detected peak was then estimated by 


/ foa \ ”pix 

P{y) = 1 “ / G{x)dx , (1) 

\'7ymax/o’y / 

where j/max is the observed peak value (minus any mea¬ 
sured continuum) in the smoothed spectrum, Uy is the 
sample standard deviation of the smoothed spectrum in 
the detection region, G[x) = e~^ 1"^ is the nor¬ 

mal distribution function, and npix is the number of 
pixels in the detection region of the spectrum. In this 
way, E(y) represents the probability of measuring one 
or more pixels with value y > j/max within the detection 
region under the null hypothesis that the pixel values 
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Fig. 1. — The bandpasses of the primary narrowband and broadband filters used to select and characterize the LAEs in this paper. Black 
curves are the four narrowband filters (from left to right: NB4325, NB4430, NB4535, NB4670, with names reflecting the central wavelength 
in Angstroms). The blue, green, and red curves show the bandpasses of the LRIS B, G, and TZ filters, respectively. 



23 24 25 26 27 

NB magnitude 


Fig. 2.— The apparent magnitude distribution of all photomet¬ 
ric (grey) and spectroscopic (black) LAEs observed in the survey. 
NB magnitude denotes the apparent magnitude of each LAE in 
the narrowband (NB) filter used to select it. The white sections of 
the bars denote the number of LAEs in each bin observed spectro¬ 
scopically for which no line was detected by our algorithm. These 
objects are excluded from the spectroscopic sample of LAEs. For 
mNB > 26.5 (hatched region), the success rate of spectroscopic 
identification of the Lya emission line dropped precipitously; for 
this reason, this paper primarily considers only those LAEs with a 
NB magnitude brighter than 26.5 (FLyQ, > 6 — 10 X erg s”^ 

cm-2). 

are normally distributed with variance cr^. Spectra with 
Piy) > 0-05 were rejected as non-detections, while those 
with P{y) < 0.05 were designated candidate emission 
lines. This threshold was chosen because it closely re¬ 
produces the samples of significant and non-significant 
spectral Lyct detections identified by eye. The median 
signal-to-noise ratio (S/N) of the spectroscopic detec¬ 
tions is 10.1. The total flux in each candidate emission 
line was then estimated both by Gaussian fitting and by 
direct integration of the unsmoothed spectrum. Lastly, 
the line was re-fit using an asymmetric line profile: 


/a (A) 


Ae A < ^ 

^g-(A-/^)/2<T2gj A > ^ 


( 2 ) 


flasym — fl'red/l^blue ■ (d) 

Here, Ubiue represents the line-width measured from 
the “blue” side of the peak (A < fi), while Ured is the 
width on the “red” side of the peak (A > /x), where /x is 


the fitted peak of the Lya profile. From these quantities, 
we define the asymmetry of the line shap e, aasym, a po¬ 


tential diagnostic of Lya escap e physics (Zheng fc Wal-| 


lace 2013 Chonis et al. 2013). Many hires were found 
to be mufti-peaked, defined by the existence of a second 
peak in the smoothed spectrum at least 2.5a above the 
continuum and within 1500 km s“^ of the primary peak, 
separated from the primary peak by an inter-peak min¬ 
imum. As these secondary peaks typically have lower 
S/N than the primary peak, they were fit using sym¬ 
metric Gaussian profiles. No more than two Gaussian 
peaks were allowed in each fit. When fitting line profiles, 
the fits were constrained such that the fit width a (both 
tTbiue and CTred for asymmetric profiles) could not fall be¬ 
low the instrument resolution {a > 1.5A). Measured line 
widths, asymmetries, an d m ulti-peaked profiles are dis¬ 
cussed in detail in Sec. 13.31 As demonstrated in that 
section, the resolution of these spectra is sufficient to 
resolve the peak width, separation, and asymmetry for 
the typical objects in our survey. However, we note that 
these spectra may not resolve spectral features on very 
small {Sv < 100 km s“^) velocity scales, such as those 
predicted by radiative transfer models and observed in 
some high-resolution spectra of bright or ne arby galaxies 
(Verharnme et al.||2008[ Martin et al.||2015|). 

The observed wavelength of the Lya emission line was 
measured via two methods: direct integration of the un¬ 
smoothed spectrum (i.e. the flux-weighted line centroid), 
and the fit value of /x (i.e. the peak of the fit asymmetric 
Gaussian profile). For multi-peaked lines, the primary 
peak is that which was detected by our first-pass line de¬ 
tection algorithm (i.e., the most significant spike in the 
smoothed spectrum). Henceforth in this paper, ^Lya.peak 
denotes the redshift of the Lya line derived from the 
asymmetric Gaussian peak, while 2;Lyo!,ave refers to the 
redshift derived f rom direct integration; these values are 
compared in Sec. |3.2| 

We used spectroscopic and photometric cuts to elim¬ 
inate non-Lya contaminants from our sample. Because 
the narrowband filters used in this survey all have central 


wavelengths Ac < 4900A, [O III] emitters are excluded by 
design. [O II] emitters (AA3726,3729) may be selected at 
z ~ 0.15 — 0.25 depending on the specific narrowband 
filter used for each field. However, the total comoving 
volume selected by these filters is ~50x larger for Lya 
at z ~ 2.6 than for [O II] at z ^ 0.2. The [O II] lumi¬ 
nosity density is also very low at z « 0.2 with respect 
to higher redshifts, and objects with [O II] equivalent 
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widths in emission greater than our o bserved-frame cut 


off (Wnhs > 727^) are extremely rare (Hogg et al. 


1998 


Ciardullo et al. 20131. Furthermore, the [U llj doublet 


should be marginally resolved in our « 100 km s“^ 
spectra, with a li ne s eparation AAqii ~ 220 km s“^. 
As shown in Sec. |3.3[ there is no detected population 
of objects that have multi-peaked emission lines with 
Aupeaks ~ 200 km s“^, so we do not expect that low- 
redshift [O II] emitters are a significant source of con¬ 
tamination to our photometric or spectroscopic samples. 

The observed spectral range for each spectrum depends 
on the location of the corresponding slit on the mask, but 
most slits were placed to cover a rest-wavelength range 
900A < a < 1500A at z ^ 2.7. As such, other contam¬ 
inants to the LAE sample were considered by searching 
for anomalous emission lines outside of the Lya detection 
window. A small number of such objects were found. All 
the spectroscopically-identified contaminants were lower- 
redshift AGN in which strong C IV AI550 (z ~ 1.9) or 
He II AI640 {z ~ 1.7) emission was detected by our nar¬ 
rowband filter. In each case, Lya, C IV, and He II emis¬ 
sion are all present with high equivalent width. Five ob¬ 
jects were found where the narrowband excess was caused 
by C IV emission, and another 5 objects were selected due 
to He II emission falling in the narrowband passband. All 
10 objects were removed from both the spectroscopic and 
photometric samples. These constraints define a sample 
of 422 spectra for 318 unique LAEs. Using the contam¬ 
ination results, we estimate the contamination fraction 
for the photometric sample to be 10/(318 -I- 10) = 3.0%. 

To facilitate comparison with samples of continuum- 
bright galaxies, we also an alyze a sample of 6 5 LBGs 
from the KBSS-MOSFIRE (Steidel et al. 2014). These 
65 galaxies are hereafter described as the KBSS LBG 
sample. Details of the object select ion and data collec¬ 
tion and reduction are described by Steidel et al. (2014) 
and references therein. The 65 KBSS LBG galaxies dis- 
cussed here are a small subsample of the entire KBSS- 
MOSFIRE catalog that meet the following requirements: 


1. rest-UV spectra acquired using the same LRIS in¬ 
strument setup (including the 600-line grism) as the 
KBSS-Lya LAEs; 


2. Lya emission spectroscopically-detected by the same 
line-detection algorithm employed for the KBSS-Lya 
LAEs (56% of the above objects); 


3. systemic redshift measurements acquired as part of 
the KBSS-MOSFIRE via rest-frame optical emission 
lines (similar to the KBSS-Lya MOSFIRE subsample 
described in Sec. 2.3); 


4. a lack of broad, high-equivalent width G IV emission 
indicative of AGN activity (97% of objects meet this 
requirement). 


The KBSS LBG sample provides a means of comparing 
the spectral properties of LAEs to objects ~10x brighter 
in the UV continuum using the same spectral resolution 
and analysis methods. Note that ~50% of all galaxies 
in the full KBSS-MOSFIRE sample have UV redshifts 
measurable only in continuum absorption lines. Such 
objects cannot be selected by narrowband searches, so 
they were omitted from the KBSS LBG sample to ensure 


appropriate comparison (criterion ^2 above). However, 
note also that the KBSS LAEs have no photometric cut 
based on Whya, as these galaxies occupy a large range of 
redshifts and do not generally have narrowband images 
probing their Lya emission. 

The line-fitting algorithm was changed slightly for the 
KBSS LBG spectra to accommodate the strong contin¬ 
uum, which often varies significantly on the redward and 
blueward sides of the Lya line. When fitting the LBGs, 
a free parameter was added to provide a linear fit to the 
continuum over the line region, but the one-peak or two- 
peak asymmetric Gaussian fits were otherwise conducted 
in an i dent ical manner for the LAE and LBG samples. 
In Sec. 4.3 we also consider a composite LBG spectrum 
from a survey of LBGs at z ^ 3 (G. Steidel et al., in 
prep.). The data collection and reduction of the objects 
in this composite spectrum are similar to the other LAEs 
and LBGs described above and will be discussed in detail 
in future work. 


QSO red s hifts w ere determined as described in Trainor 
& Steidel (2012) and have estimated uncertainties 
Cz.QSO ~ 270 km s“^. The redshift distribution of LAEs 
with respect to their nearby QSO is shown for each field 
in Fig. On small scales, the redshift distribution 
of LAEs with respect to the HLQSOs is more mean¬ 
ingful to consider as a distribution of velocities. The 
relative LAE-QSO velocities are calculated as ulae = 
c(zlae - ^QSo)/(l + -ZQSo)- The distribution of QSO- 
centric velocities for our entire spectroscopic LAE sample 
is given in Fig. Here and elsewhere in this paper we 
assume that theLAEs have systemic redshifts given by 
-^LAE = ZLya ~ 200 km s“^ to account for the typical red¬ 
shift of the Lya emission peak with respec t to systemic 
derived from our MOSFIRE data set (Sec. |2.3|). 


2.3. MOSFIRE Spectroscopic Sample 

For a small subset of these Lya-selected objects, rest- 
frame optical spectra were also obtained via the new 
Multi-Obiect Spectrometer For InfraRed Exploration 
(MOSFIRE; jMcLean et al.||2010l [20T^ on the Keck 1 
telescope. These data wer e faken over the cou rse of the 
KBSS-MOSFIRE survey (jSteidel et al. 2014). All the 
MOSFIRE data considered in this analysis are from the 
field around the HLQSO Q2343-I-12 (z = 2.573). i7-band 
and/or iL-band spectra were obtained for 43 LAEs in this 
field and reduced via the MOSFIRE-DRP; the observa¬ 
tional str ategies employed and reduction software are de¬ 
scribed in Steidel et al. (2014). Redshift fitting and one- 
dimensional spectral extraction were performed using the 
IDL program MOSPEG (A. Strom et al., in prep.). With 
OH7 slits, these spectra have R ~ 3600, ^ 35 km s“^. 

Of the 43 observed LAEs, 32 yielded robust redshifts via 
measurement of Ha (AT-band, 29 objects), H/3, and/or 
the [O III] AA4959, 5007 doublet (iL-band, 12 objects) in 
emission. Nine objects have measured redshifts in both 
bands. Four of the MOSFIRE objects do not meet the 
Lya equivalent width threshold employed here to define 
LAEs in a strict sense (WLya > 20A), but these four are 
kept in the sample because they were otherwise selected 
via the same photometric techniques used to define the 
rest of the sample, and they all display strong Lya emis¬ 
sion in their spectra. Furthermore, the other photomet¬ 
ric and spectroscopic properties of these four objects are 
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Fig. 3. — The redshift distribution of all spectroscopically-identified LAEs in each field (grey). For each field, the number of objects 
with spectra is given in the upper left. Vertical dashed lines correspond to the HLQSO redshift in each field. Solid black lines denote the 
transmission function for the narrowband filter used to select objects in each field; the normalization is such that N=1 is 10% transmission, 
the cut-off value for the automatic line-detection algorithm. Note that nearly all fields display a strong association of LAEs with the QSO 
redshift in addition to a more broadly-distributed component. 



Fig. 4. — The velocity distribution of all spectroscopically- 
identified LAEs with respect to their nearby hyperluminous QSO. 
The LAE redshifts are taken from the spectral detection of Lya, 
but are shifted by —200 km s“^ to account for the mean Lyo off¬ 
set with respect to systemic. The hatched curve is the selection 
function defined by the set of narrowband filters used to select the 
LAEs. 


typical with respect to the remainder of the MOSFIRE 
sample. 

Additional information on these spectra is given in Ta¬ 
ble and Fig. The redshift measured from the MOS¬ 
FIRE spectra, whether from the H or K band, is here¬ 
after denoted by Zneb- In the nine cases where a mea¬ 
surement was obtained in both bands, the redshift mea¬ 
surement from the higher-S/N H was used, but both 
bands generally agree quite closely: the median redshift 
difference for the two bands is 17 km s“^, roughly half 
the tJ ~ 35 km s“^ resolution of these spectra. The fit 
redshifts have typical uncertainties an ~ 5 km s“^ and 
o'K ~ 7 — 20 km s“^ for the H and K band redshifts, 
respectively, based on their Gaussian fits. 

Our pre-existing 77.-band image in the Q2343 field is 
rotated with respect to the narrowband and R-band im¬ 
ages used to select LAEs, so ~16% of the Lya-imaged 
field does not currently have 7?.-band coverage. Objects 
falling in this region are omitted when discussing aver¬ 
age 77.-band properties of the LAEs. In particular, 4 of 
the 32 MOSFIRE LAEs have no measured TZ magnitudes 
(Table but their continuum properties in other bands 
(B anoNIR where available) are consistent with those 
falling within our existing TZ coverage. The 3 -ct 77.-band 
limit is reported in Table for objects with measured TZ 
magnitudes fainter than that value. 

Velocity dispersions are measured (or constrained) for 
each MOSFIRE LAE from the Gaussian fit to the de¬ 
tected spectral line. As above, the highest-S/N line (gen¬ 
erally [O III] A5007) width is used when multiple lines 
are detected. The measured line widths (deconvolved 
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from the instrumental profile) range from cTneb = 14 ± 9 
km s“^ to 114 ± 8 km s“^ (Table § , with a median 
(o’neb) = 34 km s“^. The object with the largest fit dis¬ 
persion (Q2343-NB2807) has a two-peaked nebular line 
profile that is not well-fit by a single Gaussian; with 
the exception of this object, no LAE has a fit disper¬ 
sion larger than 73 km s“^. Ten objects have line widths 
consistent with the MOSFIRE instrument profile and are 
given as limits (dneb < 35 km s“^). The LAEs are gen¬ 
erally unresolved in our ground-based images, and fur¬ 
ther analysis of the limited HST/WFG3 imaging in t hese 
fields is fort hcoming. The similar LAE sample of |Erb| 
et al. (2014), however, have a median effective half-light 
radius a = 0.6 kpc in HST F814W photometry. Tak¬ 
ing this value as a typical size for the KBSS-Lya LAEs, 
we can estimate their d ynamical rnasses v ia the formula 
Afdyn = 5acr^/G. As in Erb et al. (2014), we assume a 
spherical geometry for these systeiris, noting that the dy¬ 
namical mass would be smaller for a disk-like geometry. 
Further discussion of the likel y ph ysical morphologies of 
these objects is given in Sec. |4(3| For the assumed ge¬ 
ometries and sizes, the KBSS-Lya LAEs would have dy¬ 
namical masses 1.4 x 10® Mq < M^yn < 8.4 x 10® 

The median mass is (M^yn) = 8.0 x 10® Mq. 

3. LYMAN-a SPECTRAL MORPHOLOGY 

The observed luminosity and spectral profile of the 
Lya line is the result of both the input distribution of 
Lya photons produced directly from recombination and 
the resonant scattering that impedes and shapes their 
escape. The net result of these processes in star-forming 
galaxies is generally seen to be the broadening of the 
Lya line and an overall shift redward with r espect to the 
systemic redshift (e.g., Shapley et al.||2003 ), an observa¬ 
tion that can be simply explained by outflowing gas that 
drives L ya photons to significant shifts in both spac e and 
velocity (Mas-Hesse et al.|2003 Steidel et al.|2010|. The 
scattering and absorption of these Lya photons likewise 
leads to severely attenuated luminosities with respect to 
early predi ctions of star-formation-ass o ciated Lya emis¬ 
sion (e.g., |Partridge & Peebles||1967t jGiavalisco et al. 
1996). The details of these effects, however, depend on 


many factors, including the dust content, column den¬ 
sity, covering fraction, and velocity distribution of the 
scattering medium. 

As the Lya line profile bears the imprint of these di¬ 
verse galaxy and gas properties, it is an important tool 
for studying galaxy evolution. However, the degenerate 
effects of these proce sses require care i n thei r interpre¬ 
tation. For insta nce, iHashimoto et al. (2013), Schenker 
et al.| (|2013|), and|bhiDuya et al.|(|2U14p show that galax¬ 
ies selected via their Lya emission have a smaller veloc¬ 
ity shift of Lya with respect to systemic compared to 
conti n uum-selected gal axies (e.g.,those of Steidel et al. 
2010). |Erb et al.|(2014| extend this study to a large sam¬ 
ple of LAEs and LBGs with systemic redshifts, demon¬ 
strating that the Lya velocity shift shows a significant 
inverse trend with WLya for the combined population of 

® With the exception of Q2343-NB2807 (for which a single Gaus¬ 
sian function is a poor fit to the nebular line profile), the largest 
estimated dynamical mass in our LAE sample is 3.7x 10® Mq . Note 
also that the lower mass bound does not include those objects with 
only an upper bound on their nebular velocity dispersion. 


galaxies selected by either technique. The physical basis 
of this relationship remains unclear, however. In partic¬ 
ular, it is not clear whether the velocity offset of Lya is 
pr edominantly driven by outflow velocity (as suggested 
by IHashi moto et al.||2013|) or H I optical depth (as sug- 
gested by Ghonis et af.||2U13 ). 

Furthermore, many LaEs do not exhibit simple Lya 
shifts, but rather show complex, multi-peaked profiles 
that may correspo nd to either internal star formation 
(Kulas et al.||201^ or exteriially-illuminated fluorescent 


processes (Kolimeier et al. 


2010). In particular, Lya 


lines that are blueshijted witn respect to systemic are pre- 
dicted by simple mode ls of emission from inflowing H I 
( Verhamme et al.|2006 ) whose low star-formation and in- 
trinsic Lya-production would likely hide them from view 
in non-fluorescing populations. The distribution of Lya 
spectral properties and their relation to other physical 
properties of the galaxies (including the Lya escape frac¬ 
tion) are thus an important window into the processes 
that drive gas into and out of galaxies throughout their 
evolution. 

3.1. Lya escape fraction 

Postage stamps of each Lya and nebular line spec¬ 
trum in the MOSFIRE sample are displayed in Fig. 
The majority of the nebular redshifts are derived from 
Ha, for which the intrinsic flux ratio is Fbya/Fna « 8.7 
under the typical assumption of case-B (i.e., ionization- 
bounded) recombination and Eg = 10^ Kj^ 

The LAE Lya line fluxes (pLya) are measured photo¬ 
metrically from their continuum-subtracted narrowband 
detections. Directly integrating the rest-UV spectrum 
yields a Lya flux ^2x smaller on average. This differ¬ 
ence may be partially due to uncertainties in the ab¬ 
solute flux calibration of the rest-UV spectra, but the 
consistency between our stacked cont inuu m spectra and 
broadband flux measurements (Sec. 14.11) suggests that 
this effect is minimal. Similarly, the photometric errors 
in the measured narrowband Lya fluxes are estimated 
to be ^15%. Rather, it appears that the differential slit 
losses between Lya and continuum photons are signifi¬ 
cant within our 172 wide slits. 

For the MOSFIRE subsample of LAEs with AT-band 
spectra, Ha fluxes (Ena; Table are estimated by a 
Gaussian fit to the emission line, slit losses are difficult 
to estimate for individual faint objects, so the measured 
values have been multiplied by a factor 1.5—2.1 (depend¬ 
ing on the mask) to account for the average slit loss mea¬ 
sured for bright point sources on the same slitmasks (A. 
Strom et ah, in prep.). Most LAEs in our sample do not 
have robust detections of both Ha and H/3, so we have 
not corrected the Ha fluxes for dust attenuation, but the 
four objects with Ha and H/3 SNR > 2 have a mean ra- 


The assumed ratio Fbya/t^Ha 


8.7 is commonly adopted in 
l(|19711l. 

4e i^T, ' 


the l iterature, often with reference to |Brocklehurst| (|197l||. How¬ 
ever, Brocklehurst (TotIJ does not explicitly calculate -^'Lya/^Ha! 
but only line ratios including states a t energy levels n > 2 (e.g., 
Hq:/H/ 3; see note in|Henry et al.|2015||. While there is some varia¬ 
tion in the calcula ted Lya/tiU and tla/lld rat ios with temperature 
and density (e.g., |Dopita &; Sutherland|2003|l, values in the range 
^LyaZ-^Ha “ ^ -- lU are typical and do not change our results 
appreciably. We adopt ^LyaZ-^Ha ~ consistency with pre- 

vious studies at h igh and low redshifts (e.g., |Hayes et al.||2010| 
[Henry et al.|20T^ . 
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TABLE 2 

MOSFIRE Observations 


Object Name 

ITLya 

(A) 

„peaka, 

Lya 


ZK^ 

peak 
'^Lya 
(km s“!) 

,,ave 

Lya 

(km s“!) 

'^neb^ 
(km s“^) 

n 

(mag) 

A,ya“ 

(10-!7) 

(10-17) 

-^Lya/ 

i^Ha 

Q2343-NB0280 

52.6 

2.5801 

2.5777 

— 

201 

364 

34±9 

26.6 

6.4±0.2 

- 

— 

Q2343-NB0308® 

58.8 

2.5690 

2.5663 

(2.5666)'^ 

226 

360 

<35 

26.4 

4.0±0.3 

1.3±0.3 

3.1 

Q2343-NB0345 

36.9 

2.5917 

- 

2.5876 

343 

414 

21±6 

25.3 

13.2±0.4 

3.8±0.3 

3.5 

Q2343-NB0405 

134.8 

2.5892 

2.5862 

(2.5860)'! 

247 

293 

33±8 

27.1 

2.9±0.2 

1.0±0.2 

2.9 

Q2343-NB0565 

313.9 

2.5661 

- 

2.5635 

215 

54 

34±12 

26.8 

5.7±0.3 

1.3±0.2 

4.4 

Q2343-NB0791 

37.4 

2.5710 

- 

2.5741 

-259 

-116 

39±9 

- 

3.6±0.3 

2.7±0.3 

1.3 

Q2343-NB1154 

54.1 

2.5773 

- 

2.5765 

65 

219 

<35 

- 

5.8±0.2 

1.0±0.3 

5.6 

Q2343-NB1174 

236.6 

2.5509 

- 

2.5478 

262 

283 

<35 

- 

12.5±0.4 

1.5±0.3 

8.6 

Q2343-NB1361 

90.6 

2.5599 

2.5587 

(2.5586)'! 

105 

50 

36±8 

27.0 

3.4±0.2 

0.3±0.1 

11.2 

Q2343-NB1386 

25.8 

2.5691 

- 

2.5647 

371 

89 

<35 

>27.3 

1.6±0.2 

1.1±0.3 

1.5 

Q2343-NB1416 

20.7 

2.5602 

2.5582 

(2.5594)'! 

169 

275 

35±14 

26.8 

1.8±0.2 

0.7±0.3 

2.4 

Q2343-NB1501 

35.2 

2.5614 

2.5593 

G.5592)'! 

180 

158 

71±12 

>27.3 

1.7±0.2 

1.2±0.4 

1.4 

Q2343-NB1518 

22.3 

2.5890 

- 

2.5860 

251 

238 

<35 

26.2 

3.3±0.2 

0.8±0.3 

4.0 

Q2343-NB1585 

303.4 

2.5669 

2.5649 

(2.5652)'! 

169 

93 

<35 

>27.3 

5.1±0.2 

0.6±0.2 

9.0 

Q2343-NB1692 

359.2 

2.5625 

- 

2.5602 

194 

85 

14±9 

26.8 

5.6±0.2 

1.1±0.3 

5.0 

Q2343-NB1783 

35.0 

2.5784 

2.5767 

(2.5766)'! 

143 

167 

51±3 

25.7 

7.8±0.2 

2.0±0.2 

3.8 

Q2343-NB1789 

88.2 

2.5471 

- 

2.5448 

196 

59 

45±16 

>27.3 

2.3±0.2 

1.4±0.3 

1.6 

Q2343-NB1806 

19.3 

2.5982 

- 

2.5954 

237 

204 

73±19 

26.7 

1.2±0.1 

1.4±0.3 

0.8 

Q2343-NB1828 

26.7 

2.5753 

- 

2.5724 

247 

216 

38±15 

26.9 

2.7±0.3 

1.0±0.2 

2.6 

Q2343-NB1829 

32.1 

2.5782 

- 

2.5754 

233 

271 

41±8 

25.9 

5.7±0.2 

2.6±0.3 

2.2 

Q2343-NB1860 

23.6 

2.5740 

- 

2.5741 

-11 

237 

<35 

25.7 

1.2±0.2 

1.0±0.3 

1.2 

Q2343-NB2211® 

66.2 

2.5776 

2.5760 

(2.5758)'! 

135 

206 

30±13 

26.9 

4.3±0.3 

1.1±0.2 

4.0 

Q2343-NB2571 

29.4 

2.5820 

- 

2.5850 

-249 

-503 

38±26 

>27.3 

1.4±0.3 

0.8±0.2 

1.8 

Q2343-NB2785 

14.2 

2.5723 

- 

2.5785 

-517 

-367 

<35 

>27.3 

1.3±0.1 

0.4±0.1 

3.2 

Q2343-NB2807'= 

24.6 

2.5479 

2.5446 

- 

279 

184 

110±8 

23.8 

23.4±0.3 

- 

- 

Q2343-NB2816 

58.5 

2.5782 

- 

2.5742 

340 

243 

57±23 

>27.3 

2.0±0.1 

1.5±0.4 

1.4 

Q2343-NB2821 

32.4 

2.5800 

- 

2.5779 

173 

69 

24±23 

- 

1.3±0.2 

1.3±0.4 

1.0 

Q2343-NB2834 

7.2 

2.5700 

- 

2.5683 

142 

-69 

<35 

>27.3 

0.9±0.1 

0.5±0.2 

1.7 

Q2343-NB3061 

57.6 

2.5788 

- 

2.5764 

203 

270 

72±24 

26.5 

2.2±0.2 

2.0±0.5 

1.1 

Q2343-NB3170 

90.0 

2.5488 

2.5475 

- 

113 

152 

15±5 

>27.3 

2.6±0.1 

- 

- 

Q2343-NB323lf 

110.4 

2.5737 

2.5711 

(2.5695)'! 

218 

151 

43±12 

25.9 

8.4±0.2 

1.0±0.2 

8.7 

Q2343-NB3292 

13.4 

2.5700 

- 

2.5631 

579 

153 

<35 

>27.3 

1.3±0.2 

0.6±0.2 

2.2 


^ denote the redshifts measured from fitting the Lya emission line, H band emission lines (primarily [O III] A5007), 

and K band emission lines (Ha exclusively). 

^ cTneb i® velocity width of the nebular emission line used to define Both the line widths and their associated uncertainties 

reflect the deconvolution of the = 35 km s~^ instrumental profile. 

^ Line fluxes (in cgs units) are determined from a combination of spectroscopic and photometric measurements as described in Sec. |3.2[ 
^ In the nine cases where both the H and K band spectra yielded redshifts, we adopted the higher-S/N H band redshifts to define -^neb- 
® The nebular redshifts of Q2343-NB0308, Q2343-NB2211, and Q2343-NB2807 were measured from the [O III] A4959A line because the 
A5007A line did not fall on the detector. 

^ The nebular redshift of Q2343-NB3231 was measured from the H/3 line (and cross-checked with Ha) because the [O Hi] A5007A and 
A 4959 A lines did not fall on the detector. 


tio (Fhq/-Fh/ 3 ) = 3.3zb 1.1, consistent with the ca nonical et al. 


(non-attenua ted) case-B value Fua/Fiijs = 2.86 (Brock- 
lehurstlpml ). 


i'he measured line ratios in Fi^[^ vary consider¬ 
ably {Fi^ya/Fna 1 — 11, Table H with respect to 
that predicted by case-B recombination. If the intrinsic 
Fhya/Fiia Aux ratio is assumed to be 8.7 (and the Ha 
flux is not significantly affected by extinction), then the 
median ratio {Fi^ya/Fua) =2.7 suggests a typical Lya 
escape fraction (/esc.Lya) = « 30%. 

This Fi,ya/F}ia ratio is consistent with that observed 
in the 4 objects wit h Ha line detecti ons in the similar 
faint-LAE sample of Erb et al. (2014). If the Ha emis¬ 
sion is attenuated by interstellar dust, then the intrin¬ 
sic values of /esc,Lya may be lower, but the measured 
values are considerably higher than previous estimates 
of .fesc.T.y g ~ 5% for typical high-redshift star-forming 
galaxies ( Hayes et al.|20io| Steidel et al.|20lT Ciardullo 


2014). Such high Lya escape fractions are sugges- 


tive of the high LyC escape fractions observed in sam ples 
of LAEs (Nestor et al. 2013 Mostardi et al. 2013) and 
may indicate low H 1 optical depths tor photons escaping 
from star-forming regions in these galaxies. Given the 
proximity of these LAEs to hyperluminous QSOs, how¬ 
ever, it may be that /esc,Lya is elevated by the “fluores¬ 
cent” boosting of the total Lya flux from these objects. 

Eor comparison, the KBSS LBGs have a median flux 
ratio (ELya/EHa) ~ 0.5, or a corresponding escape frac¬ 
tion /esc,Lya ~ 6%- The LBG Ha slit losses were cor¬ 
rected in a similar manner to the LAEs, and further anal¬ 
ysis of the rest-frame optical spectra of the KBSS LBGs 
will be presented in A. Strom et al., (in prep.). The LBGs 
lack narrowband Lya photometry, so the Lya fluxes were 
measured by direct integration of the rest-UV spectrum 
(after flux-calibratin g the spectrurn based on the broad¬ 
band G-band flux). jSteidel et al. (2011) measure Lya 
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aperture corrections for LBG spectra showing Lya in net 
emission or absorption in comparison to LAEs, finding 
that LBGs with detected Lya emission show similar dif¬ 
ferential slit losses to LAEs in our l'/2 slits. In this case 
the Lyof luminosities of the KBSS LBGs may be under¬ 
estimated by a factor of ^ 2 x, and the escape fractions 
may be /esc,Lya ~ 12% in a larger aperture. This Lya es¬ 
cape fraction is similar to t hat derived fo r Lyg- emitting 
LBGs (/esc.Lya = 14.4%) by 
appears that our sample o: 


visible. This shift is consistent with that found for the 
few samples o f LAEs with systemic redshifts in the lit- 


Steidel et al. (2011). It thus 
Lya-emitting, continuum- 
color selected KBSS LBGs have a broadly similar Lya 
escape fraction with respect to previously studied sam¬ 
ples of LBGs, and significantly smaller values than the 
KBSS-Lya LAEs. 

Given the small number of LAEs with Ha detections, 
it is difficult to investigate trends within the MOSFIRE 
LAE sample. In particular, we find no association be¬ 
tween /esc and the Lya offset, peak multiplicity, or peak 
separ ation (uT. vg, /mult, and Upeaks) see discussion in 
Secs. 3.2-3.31. However, we find a significant trend be¬ 
tween j^cLya and Wkya, which are positively correlated 
with p < 4 X 10“^ for a Pearson correlation test. This re¬ 
sult suggests that WLya is primarily an indicator of Lya 
escape, not Lya production among these faint galaxies. 
In fact, the relation between Ha luminosity (a standard 
indicator of star-formation rate) and Wj^ya is weakly neg¬ 
ative in the sample. Given that a higher Wkya at a given 
observed Lya luminosity corresponds to a fainter UV 
continuum luminosity (another proxy for star-formation 
rate), we suggest that Wkya is likely anti-correlated with 
star-formation rate among these faint LAEs. This ef¬ 
fect is similar to the anti-correlation observed between 


in several previous studies (e.g 

, Ando et al.||2006l 

Ouchi 

et al. 

120081 IStark et al.||2010l 

Kornei et ai.| 2010 

|Atek 

et al. 

|2014|). We also find a tentative, negative associa- 


and/or [O III] emission). Because (Tneb traces the stellar 
velocity dispersion and thus the dynamical mass of the 
galaxy, this may suggest that higher WLya is associated 
with lower dynamical masses. However, further study of 
the stellar populations and star-formation rates of the 
LAEs (including our ongoing photometry and rest-frame 
optical spectroscopy of these objects) will be necessary 
to confirm the physical basis of these trends. 

3.2. Velocity shift with respect to systemic redshift 

The offset of the Lya line with respect to the systemic 
redshift was measured for each of the 32 LAEs in the 
MOSEIRE sample using the Zuy g andzneb values esti¬ 
mated as described in Sec. |2.2| and |2.3[ The velocity 
offset was then inferred from the measured redshifts as 


^Lvci ^neb 
VLya = -C . 


^neb 


( 4 ) 


The velocity shift of the Lya line with respect to the 
systemic redshift is clear from the panels in Fig. 5 The 
majority (28/32; 88%) of spectra display the redshifted 
Lya profile typical of star-forming galaxies, but there 
are several objects with minimal velocity shift, or even 
a dominant peak blueward of the systemic redshift. The 
distribution of measured velocity shifts (Eq. 13 is given 
in Fig.[^ The typical Lya shift of -1-200 km s“^is clearly 


erature (e.g., McLinden 

et al. 20111 iHashimoto et al. 

2013| Ghonis et al. 

2013 

Shibuya et al. 

2014), and sig- 


nihcantly less than that previously measured in samples 


LBGs (300—400 km s Steidel et al. 2010). However 


we demonstrate below that the detailed velocity distribu¬ 
tion depends on the method by which the Lya emission 
redshift is assigned. 

The two panels of Fig. show the effect of defin¬ 
ing ^Lya by either the fitted peak (^Lya.peak) or the 
flux- weig hted centroid of the line profile (zLya.ave) see 
Sec. |2.2[ ). Both methods yield qualitatively similar dis¬ 
tributions, with median values (wLya,peak) = 198 km s“^ 
and (uLya.ave) = 176 kms“k However, zuya.peak displays 
a tighter correlation with the redshift derived from nebu¬ 
lar emission lines than ZLya.ave- The median absolute de¬ 
viations (a scale estimator that is insensitive to outliers) 
of the two distributions of VLya are MAD(uLya,peak) = 54 
km and MAD(uLya,ave) = 93 km The similar¬ 
ities between the two distributions suggest that either 
quantity is a reasonable proxy for the systemic redshift 
after correcting for the typical 200 km s“^ offset, and 
even the individual velocities computed by both measure¬ 
ments are consistent within 100 km s“^ for the majority 
of objects (Table [^. Because ZLya,peak tracks Zneb more 
tightly, however, we adopt Zsys,Lya = ^Lya.peak - 200 as 
the best estimate of the systemic redshift in the absence 
of nebular emission line measurements. 

The distribution of Lya line offsets from the compar¬ 
ison KBSS LBG sample is also displayed in Fig. and 
there the differences between the two Lya redshift esti¬ 
mators are more apparent. While the distribution of fit 
peak velocities ULya,peak is again sharply peaked around 
VLya ~ 200 km s“^, the centroid velocities WLya.ave 
are more extended, particularly toward higher (more 
strongly redshifted) velocities. The median velocity off¬ 
sets for the LBG sample are (uLya,peak) = 205 km s“^ 
and (uLya,ave) = 306 km s“^. Note that our measurement 
of ULya,peak in particular suggests a significantly smaller 
Lya offset than previous studies of LBGs. Qualitatively, 
the difference between the two distributions is due to the 
fact that the LBG Lya emission lines have extended red 
tails (often combined with continuum flux levels that are 
higher on the red side of Lya than the blue side) that 
draw the velocity centroid to higher/redder values. Fit¬ 
ting the Lya emission line with a symmetric Gaussian 
model will likewise cause the fit emission-line center to 
be redder than the Lya emission peak. While the LAE 
emission lines tend to be asymmetric as well, they typi¬ 
cally have more symmetric profiles and/or stronger blue 
secondary peaks that compensate for the extended red 
wing. The LAEs also have less continuum emission that 
is more symmetric about the wavelength of Lya than 
typical LBG s (e.g., the observed ratio of /1125 to / 1325 ; 
see Sec. 4.3). The distributions of emission line widths, 


asymmetry, and peak multiplicity for each sample are 
discussed in Sec. 13.31 below. 

Despite their differing line morphologies and photo¬ 
metric properties, it is significant that both the LAE 
and LBG samples have Lya emission line peaks that are 
well-described by a uniform offset of ^200 km s“^ red- 
ward of the systemic redshifts. While the measured ve- 
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resolution R ~ 1300. Nebular spectra are for the Ha line (K band, red) or H/3 or [O III] {H band, purple), whichever was used to estimate 
the systemic redshift (see Tableand notes below). MOSFIRE spectra have been smoothed by a three-pixel boxcar to reduce noise. Grey 
regions denote the itlcr uncertainties on the MOSFIRE spectra. Both spectra are shifted to the redshift frame estimated from ^neb- 'll*® 
fluxes are in XFx units to facilitate comparison of their integrated line luminosities. (Figure is continued below) 


locity of the Lya peak may be strongly dependent on 
spectral resolution, its constancy across a broad range 
of galaxy luminosities suggests that it is less sensitive to 
galaxy properties than the Lyct emission centroid. For a 
sample of galaxies observed with the same spectral reso¬ 
lution, therefore, the Lya emission peak appears to be a 
robust redshift indicator across a diverse range of galaxy 
types (with the caveat that ~50% of z ^ 2.4 LBGs hav e 
no net Lya emission in their spectra; Steidel et al.|2011|). 
We use this redshift calibration to create composite LAFl 
spectra from our large sample of LAEs without nebular 
redshift measurements in Sec. U) 

In addition to the distribution of average or peak Lya 
line velocities, it is interesting to compare the full veloc¬ 
ity distributions of Lya and nebular emission. Fig. 0d is- 


plays the stacked Lya and nebular line profiles of the 32 
MOSFIRE LAEs. Of particular note is the breadth of the 
Lya line with respect to the nebular line profile; assum¬ 
ing that the Lya photons are generated by recombina¬ 
tion processes in the same H II regions that generate the 
nebular emission, the top panel of Fig. [^clearly displays 
the diffusion of the Lya photons in velocity as they reso¬ 
nantly scatter through the surrounding H I gas. Fitting 
a Gaussian function to the stack of nebular lines yields a 
velocity width of Cneb = 35±3 km s“^. The stacked Lya 
profile has a width of cr^ya = 115 ± 8 km s“^ for the pri¬ 
mary (red) peak. These figures reflect the deconvolution 
of the instrumental resolution of MOSFIRE (crinst = 35 
km s“^) and LRIS (cinst = 100 km s“^) in their observed 
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Fig. 5 (CONT.).— Spectra for the remainder of the MOSFIRE LAE sample. 




modesn The peak of the stacked Lya spectrum is again 
quite close to +200 km s“^, but there is a distinct com¬ 
ponent blueward of the systemic redshift as well. After 
subtracting the faint continuum, the fraction of Lya line 
flux emitted with n < 0 is 16%. For comparison, we 
show the effective profile that would result from stacking 
all the Lya lines at the redshifts derived from the Lya 
line peak and then making a +200 km s“^ shift in the 
bottom panel of Fig. The resulting profile is narrower 
(o'Lyct = 85 + 8 km s^), with only 11% of the line flux 
emitted at z; < 0, roughly 2/3 the flux measured using 
the true systemic redshifts. 

The four MOSFIRE LAEs with VLya < 0 compose 
an extremely interesting, though small, group of ob¬ 
jects. Further analysis of this sample is ongoing, includ- 

® Note that the fits to the composite line profiles yield similar 
(though not identical) values to the averages of the individually-fit 
lines (Table 1^. 


ing spectroscopy and photometry to constrain the stel¬ 
lar contribution to the LAE luminosity, enrichment, and 
dynamics. While only three of these blueshifted LAEs 
currently have TZ photometry, their average flux in that 
band is ~l/2x that of the LAEs with redshifted Lya 
lines, tentatively suggesting that blueshifted Lya emit¬ 
ters have physical properties that are distinct from typ¬ 
ical LAEs. Other properties associated with redshifted 
or blueshifted Lya emission are discussed in the context 
of multi-peaked Lya emission below. 


3.3. Multi-peaked and asymmetric Lya profiles 

As discussed above, star-forming galaxies are often as¬ 
sociated with multi-peaked Lya emission. In a system¬ 
atic study of Lya emis sion among ^1500 star-forming 
galaxies at 2 ~ 2 — 3, Kulas et al. (2012) found that 
^30% are multi-peaked, with some dependence on the 
spectral resolution and S/N of the observations. While 
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Fig. 6. — Distribution of Lya velocity offsets with respect to 
systemic for LAEs (grey) and KBSS LBGs (orange) with nebu¬ 
lar redshifts measured with MOSFIRE. In each panel, the vertical 
dashed line denotes ULyg, = 0. Top: The distributions when zpya 
is defined by the peak of the fit line profile (zpyo! peak)- Most of 
the galaxies have offsets in a tight range around = 200 km 

s“^, but a few LAEs have Lyo lines blueshifted with respect to 
the systemic redshift (upy^j < 0). Bottom: The distributions when 
^Lya defined by the line centroid (zpya ave)- Most LAEs have 
^Lya ~ 200 km s“^, but with a broader distribution compared to 
the zpya = -^peak distribution. The LBG centroid velocities oc¬ 
cupy a broader, more highly-redshifted range than either the LAEs 
or the corresponding LBG Lya peaks. 


most of these galaxies were observed at lower resolution, 
a subsample of their spectra were obtained with the same 
600-line grism of Keck/LRIS used for our LAE sample 
and have similar integration times (~1.5 hours). This 
subsample included 44 multi-peaked spectra, a multiplic- 
it y rate of 27%. 

(|2012a|) studied the Lya peak mor- 
;AEs at z ^ 3 at similar resolution 


Yamada et al. 
phologies of 91 
{R ^ 1700), finding a ~50% multiplicity rate. This het- 
erogeneous s ample contains 12 Lya-blobs from |Matsuda| 


et al. (2004), as well as compact, faint LAEs similar to 
those of the KBSS-Lya sample (although with a limit¬ 


ing Ly g flux ~2x brighter, according to [Yamada et al. 


2012b). 

Tor the KBSS-Lya sample of this paper, multi-peaked 
systems were identified as described in Sec. |2.2| Of 


the 318 unique spectroscopically-identihed LAEs, 129 are 
found to be multi-peaked, for a multiplicity fraction of 
40%. The criteria for id entifying a seconda ry peak were 
not identical to that of Kulas et al. (2012); in particu¬ 
lar, the significance threshofd was not dependent on the 
magnitude or S/N of the primary peak, and no mini¬ 
mum peak separation was enforced. However, we found 
that adjusting the multiplicity criteria to match those of 


Kulas et ^ as closely as possible did not significantly af- 
fect the overall rate of selection nor the spectral/physical 
properties of the selected sample. 

Previous studies of Lya kinematic multiplicity (among 
samples with or without systemic redshifts) typically 
group these lines into blue-peak dominant or red-peak 
dominant lines, motivated by the association of blue 




Fig. 7.— Stacked spectral profiles for the 32 LAEs with sys¬ 
temic (nebular) redshift measurements. Top: Comparision of the 
stacked nebular (red) and Lya (black) line profiles for these LAEs. 
The nebular stack consists of all measured nebular lines (Ha, H/3, 
and/or [O III] AA4959,5007) for each LAE. The lines are stacked 
according to their corresponding nebular redshifts, and the effect 
of resonant scattering on escaping Lya photons is clearly visible, 
as is the typical Lya velocity offset of ~ +200 km s~^ relative to 
systemic. Bottom: Comparison of the average Lya profiles when 
spectra are stacked according to their nebular, systemic redshift 
(black, as above) and a redshift derived directly from the Lya line 
peak (blue, shifted +30 km s~^ for clarity), wherein the “systemic” 
redshift is estimated as Ugyg Ly^ = 'i^Lyo, — 200 km s“^. Stacking 
via the Lya redshift distorts the Lya profile and diminishes the 
measured flux blueward of the systemic redshift. 
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Fig. 8. — A sample of the diverse Lya spectral morphologies represented within our spectroscopic sample of LAEs. Red vertical lines 
denote the redshift fit to the primary Lya peak by our algorithm, while blue vertical lines denote the fit position of the secondary peak 
(when detected with significance). Redshifts are assigned to be —200 km s“^ with respect to the primary Lya peak. The top row (panels 
a—d) includes examples of common morphological types: single-peaked lines with red (a) or blue (b) asymmetric tails, or double-peaked 
lines with dominant red peaks (c, d). The bottom row (panels e—f) shows examples of the less common blue-dominant double-peaked 
lines. Despite their similar peak ratios, panels (e & /) have extended blue tails, whereas panels (g & h) have extended red tails. 


(red) peak dominance with infl owing (outflowing;) gas in 
radiation transfer models (e.g ., Zheng &: Miralda-Escude 


2002| iVerhamme et al. 2006 Dijkstra et al. 2006” Ku- 
las et aT ( |2012 ) found that 67% of their 239 multi- 
peaked objects (the majority of which had lower spec¬ 
tral resolution than the KBSS-Lya sample) have dom¬ 
inant red peaks. A selection of examples of the di¬ 
verse Lya spectral morphologies represented among the 
KBSS-Lya LAEs is displayed in Fig. and the dis¬ 
tribution of red/blue peak dominance in our sample is 
displayed in Fig. In the top panel of Fig. we 
show the distribution of inter-peak velocities, defined as 

Axipeaks — c(2^peak,pri ^peak,sec)/(l “t (^peak,pri)- ^peak,pri 

is the redshift of the first peak ide ntifie d by the line de¬ 
tection algorithm discussed in Sec. 2.2 (i.e., the peak of 
the asymmetric Gaussian profile fit to the highest peak 
in the smoothed spectrum), and Zpeakjsec is the redshift 
corresponding to the mean of the second fit Gaussian 
component. For this definition, Aupeaks > 0 corresponds 
to red-dominant systems, such that the primary peak is 
redward of the inter-peak trough. 

Of the 129 multi-peak LAEs, 96 (74%) have domi¬ 
nant red peaks (Aiipeaks > 0), while 33 (26%) are blue 
dominant (Aupeaks < 0; Table The red-dominant 
LAEs have a quite narrow distribution of peak separa¬ 
tions, with a median value (Aiipeaks) = 560 km s“^ and 
a standard deviation cr(Awpeaks) = 220 km s“^. Con¬ 
versely, the blue-dominant LAEs have a larger velocity 
shift over a broader range: (Aupeaks) = —660 km s“^ 
and (T(Atipeaks) = 300 km s“^. A similar contrast can be 
seen in the flux distribution within these multi-peaked 
LAEs (Fig. 1^ bottom panel). The fraction of total Lya 
flux in the red peak, Fi.ed/Ttot 7 is defined as the integral 
of the redder of the two Gaussian components divided 
by the integral of the sum of the two components. There 
are 101 LAEs (78%) with Fred/Ttot > 0.5, an alterna- 




Fig. 9. — Top: Distribution of peak separations among the 129 
LAEs with identified multi-peaked Lya emission lines (40% of the 
318 LAEs with Lya spectra). The majority (96/129) have domi¬ 
nant red peaks (An > 0, red), but a substantial minority (33/129) 
have dominant blue peaks (Au < 0, blue). The dashed grey line 
denotes An = 0. The red-peak-dominant systems have a smaller 
median peak separation (560 km s“^ vs. 660 km s~^) than the 
blue-peak-dominant systems, as well as a smaller scatter about 
this value (220 km s“^ vs. 300 km s“^). Bottom: The distribu¬ 
tion (stacked histogram) of flux in the red peak as a fraction 

of total Lya flux (Etot) for fir® same 129 multi-peaked LAEs, with 
colors matched to the top panel. Red-peak-dominant Lya lines 
(red. An > 0 in the top panel) generally have Ej-g^^/Etot > 0-6, 
and the median flux ratio among all systems is Ej-g^/Etot = 0.68 
(dotted black line). 
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Fig. 10.— The redshift distribution of multi-peaked hya lines 
with respect to their corresponding HLQSOs. Color-codi ng of the 
stacked histogram matches Fig. ^(top panel) and Fig. |11[ The 
dashed grey line corresponds toTjAEs at the HLQSO redshift. 
The dotted black line corresponds a distance 10 Mlyr behind the 
HLQSO, the upper bound on the region that could emit detectable 
QSO-gene rated fluorescent emissio n for a 20 Myr QSO lifetime (as 
argued in [Trainor fc Steidel|[^13[ |. Here, the LAE systemic red- 
shifts are assumed to be —200 km s“^ from the red Lya peak. 
Under this assumption, the blue-peaked LAEs are highly biased 
toward the rear of the nearby QSOs. Black hatched bars denote 
the systemic redshifts measured for the 4 LAEs with blueshifted 
Lya peaks peak < 0 km s“^), which appear to occupy a 

similar physical volume as the blue-dominant multipeaked LAEs. 


tive definition of red-peak dominance. The sets of LAEs 
with A^eaks > 0 and Ei-ed/Ttot > 0.5 overlap almost 
entirely^ but we use the prior definition (based on the 
significance o f peaks in the smoo thed spectrum) for con¬ 
sistency with Kulas et al. (20121 and other surveys that 
assume different model prohles. 

Due to the presence of nearby hyperluminous QSOs in 
these LAE fields, it is important to discern what role the 
external illumination of the QSO (i.e., “fluorescence”) 
may play in shaping the observed spectral morphology 
of the LAEs. Clearly, the fluorescing LAEs must inhabit 
the region of the environment illuminated by the QSO in 
order for the ionizing field to have an observable effect. 
As argued in [Trainor & Steidel| (|2013|) , the combination 
of geometry and the finite QSO lifetime (measured to 
be <20 Myr) cause this fluorescence to be observable 
only when the LAE is in the foreground of the QSO, or 
less than ~10^ lightyears (3.2 proper Mpc) distant in the 
background. 

The redshift distribution relative to the QSOs of all 
the multi-peaked LAEs is shown in Fig. 10 Care must 
be taken in comparing the redshift distributions of red- 
and blue-dominant LAEs based on the redshifts of their 
Lya lines. Fig. shows the redshift distribution when 
the systemic redsmft of each LAE is taken to be 200 km 
s“^ blueward of the red Lya peak, regardless of whether 
the redder peak is primary or secondary. The redshift 
distribution of red-dominant Lya profiles is fairly sym¬ 
metric, with a slight bias toward the QSO foreground. 
The blue-dominant LAEs, however, show a strong as¬ 
sociation with the region in the immediate background 
of the QSO. A two-sample KS test finds only a small 
{p « 0.03) likelihood of both samples being drawn from 
the same distribution. 

Given the uncertainty in assigning redshifts to 


® The small number of objects with Fred/^tot > 0-5 and 
Aupgaks < 0 (or vice versa) generally have a secondary peak that 
was fit with a broad Gaussian profile, thereby including more total 
flux despite having a lower S/N. 


Fig. 11. — Lyo line asymmetries for LAEs with multi-peaked Lya 
detections. The reported asymmetry corresponds to the primary 
peak of the multi-peaked line. Red bars give the asymmetry dis¬ 
tribution of red-peak-dominant systems (Aupgaks > 0)i while blue 
bars correspond to blue-peak-dominant systems (Aupg^j^g < 0) ac¬ 
cording to the top panel of Fig.[^ The blue-white dashed line is the 
median asymmetry of the Aupg^j^g < 0 sample and the red-white 
dashed line is that of the Aupg^j-g > 0 sample. For comparison, 
the black dashed line gives the median asymmetry of the sample 
of spectral Lyo lines without multiple detected peaks. 


those sources without accurate systemic redshifts from 
optically-thin nebular emission lines, it is unclear 
whether the blue-dominant LAEs are localized within the 
QSO held, or whether they differ from the red-dominant 
distribution at all. As noted above, blueshifted Lya 
emission is a possible observational signature of infalling 
gas illuminated by a central source. Given the proximity 
of the HLQSOs, however, it is also possible that the de¬ 
tected emission has a signihcant huorescent contribution. 
In this case, the blueshifted Lya emission could con¬ 
ceivably correspond to outflowing gas illuminated exter¬ 
nally by the nearby QSO, a situation that could explain 
the association of blue-dominant LAEs with the back¬ 
ground of the QSO, where the geometry would naturally 
predict blueshifted emission from the QSO-facing side 
of a gaseous outhow. Only three of the blue-dominant 
LAEs currently have MOSFIRE spectra, but their sys¬ 
temic redshifts place them 0.4—1.2 pMpc in the back¬ 
ground of the QSO (consistent with the Lya-only mea¬ 
surements). Similarly, the 4 LAEs with MOSFIRE red¬ 
shifts and WLya.peak < 0 km s“^ (that is, objects with 
blueshifted Lya with respect to systemic) lie in the im¬ 
mediate back grou nd of the QSO (0.4—4 pMpc; hatched 
region in Fig. 101, suggesting that external illumination 
may play a role in their spectral morphology as well. 

Detailed predictions for fluorescent emission from gas 
with large bulk velocities and a realistic spatial distribu- 
tion are unclear (although see work by Kollmeier et al.| 


20101, and it is difficult to draw robust conclusions with- 
out access to the systemic redshifts of these LAEs. Given 
the diversity of line morphologies and asymmetries even 
within the set of blue-dominant LAEs (Fig. [^, it is un¬ 
clear whether these atypical profiles are generated by a 
single mechanism. Further study of these “blue” LAEs, 
including additional rest-frame optical spectroscopy, is 
required to understand the physics of their Lya emission 
and transmission. 

Lastly, in the absence of systemic redshift measure¬ 
ments and/or multiple emis sion peaks, the Lya line 
asymmetry (aasym; Sec. 2.2) may al so trace gas kine¬ 
matic s among these LAEs (see, e.g., [Zheng & Wallacej 

;en not evident in 


2013). However, line asymmetry is of 
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Fig. 12.— Distribution of fit widths (Eq.l^ for the primary Lya 
peak among the 318 KBSS-Lya LAEs (gre^ and the 65 LBGs in 
the KBSS LBG (orange) sample lines. The LRIS instrumental res¬ 
olution of = 100 km s~^ has been subtracted in quadrature; 

fit widths less than 100 km s~^ are plotted as zero. The left panel 
shows the number distribution of object line widths, while the right 
panel presents the same data normalized as a probability distribu¬ 
tion for ease of comparison. The median width for each sample is 
marked with a dashed line. 
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Fig. 13. — Distribution of Lya peak separations among the 129 
LAEs in Fig. (grey) and the 29 LBGs in the KBSS LBG (orange) 
sample with multi-peaked Lya emission lines. Given the smaller 
size of the LBG sample, positive and negative separations of the 
same size are grouped together (unlike Fig. 1^. The left panel 
shows the number distribution of of objects, A^ile the right panel 
presents the same data normalized as a probability distribution 
for ease of comparison. The median separation for each sample is 
marked with a dashed line. 

medium-resolution, low-S/N spectra. Not only will the 
measured line asymmetry be systematically diminished 
by the smoothing effect of the convolved line spread func¬ 
tion, but noise peaks and unresolved emission line fea¬ 
tures may be fit by an extended tail by our automatic 
fitting algorithm. 

Despite these issues, some trends can be seen in the 
relative distributions of Lya line asymmetries. Fig. 
shows the distribution of aasym for blue-dominant 
red-dominant multi-peaked LAEs. While both subsam¬ 
ples have broad distributions of aasym because of the ef¬ 
fects described above, there is a clear tendency of the 
blue-dominant Lya lines to have lower values of aasym 
(that is, an extended blue wing in the primary peak), 


while red-dominant LAEs typically have aa 


> 1 


(extended red wings). The two distributions have a 
Kolmogorov-Smirnov probability of being drawn from 
the same underlying distributio n p < 4 x 10~^. Thi s 
result agrees with the findings of|Yamada et al. (2012a), 
who found a similar relation (albeit with lower signin- 
cance) in their smaller sample. 


TABLE 3 

LAE AND LBG Emission-Line Properties 



LAEs 

MOSFIRE 

LAEs=" 

KBSS 

LBGs 

-^obj 

318 

32 

65 

(l^Lya I'" (10-1^ cgs) 

4.5±0.8 

4.6±0.8 

7.3±1.6 

{(r-Lya) (km s 1) 

144±3 

138±8 

195±7 

'^blue'" 

— 

4 

4 

/blue" 

- 

13% 

6 % 


129 

13 

29 

/mult‘s 

41% 

41% 

45% 

'^mult,blue 

33 

3 

6 

f mult,blue 

10 % 

9% 

9% 

(I'^^peaksl) (km s"!) 

615±14 

608±49 

716±35 

(%a) (10-1'^ cgs) 

- 

1.3±0.1 

7.2±0.7 

(irneb} (km s 1) 

- 

34±4 

84±4 


Parameters for the subset of LAEs with MOSFIRE neb¬ 
ular (systemic) redshifts 

Averages are given as mean values, with uncertainties 
denoting the standard error of the mean. 

Number (or fraction) of objects with a blueshifted Lya 
line (uLyo, < 0) 

Number (or fraction) of objects with multiple detected 
Lya peaks 

® Number (or fraction) with multiple detected Lyo peaks 
that have a lower-velocity (“blue”) dominant peak 


We compare the distribution of Lya morphologies be¬ 
tween the KBSS-Lya LAEs and the KBSS LBGs in 


Figs 

29/65 


as well as in TableIn the LBG sample 
ave a detected secondary peak. Of these 
objects, 23/29 (79%) are red-peak dominant. Note that 
these rates are 40% and 74%, respectively, among the 
LAEs (Table [^. As described above, measured rates of 
Lya peak multiplicity are dependent on both the resolu¬ 
tion and S/N of the spectra. The KBSS LBGs considered 
in this paper were observed using the same instrument 
setup as the KBSS-Lya LAEs and were anal yzed using 
the same line-detection software (see Sec. 2.2); however, 
the Wpical LBG Lya fluxes are significantly brighter (Ta¬ 
ble 1^ and the spectral S/N is correspondingly higher. 
Secondary peaks of low relative significance are there¬ 
fore more likely to be detected in the LBG sample than 
in our LAE sample. In particular, we find that the ra¬ 
tio of primary-peak flux to secondary-peak flux is larger 
among the LBGs than the LAEs, which may be caused 
by missing similarly faint secondary Lya peaks in the 
LAE spectra. However, radiative transfer models of Lya 
escape also predict that the primary/s econdary flux ra¬ 


tio in creases with outflow velocity (e.g., Verhamme et al. 


20061. In either case, the frequency of Lya peak multi¬ 


plicity and blue-dominant Lya peaks are broadly consis¬ 
tent between our LBG and LAE samples. 

However, the velocity distribution of Lya flux varies 


significantly between the LAE and LBG samples. Fig. 12 
displays the distribution of fit Lya emission line widths 
for both samples of objects. The measured line width is 
that of the primary Lya peak, and is given by the average 
of the red and blue scale parameters of the asymmetric 
Gaussian fit (Eq. : 


^Ly O', primary 


*^red “1“ ^blue 
' ~2 


( 5 ) 


3.4. LBG Lya line morphologies 


The LAE sample has a median line width of 166 km 
s“^, with 33% of spectra having a primary-peak line 
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width less than 150 km s“^. The LBG sample has a 
median line width of 210 km s“^, with only 3% (two 
spectra) having a primary-peak line width less than 150 
km s“^. The LRIS spectral resolution crinst = 100 km s“^ 
(which is the same for both the LAE and LBG samples) 
has been subtracted in quadrature. Note that nearly 
all the measured lines are well-resolved with respect to 
the resolution limit. A Kolmogorov-Smirnov test gives a 
probability p « 10“^*^ of both samples being drawn from 
the same distribution. 

A similar contrast is seen in the distribution of peak 
separations among the LAEs and LBGs with multiple 
detected emission peaks (Eig. 131. Among the 129 LAEs 
with multiple detected peaks, the median peak separa¬ 
tion is (|Aupeaks|) = 565 km s“^, while the 29 LBGs with 
multiple peaks have (|A?;peaks|) = 780 km s“^. Note that 
unlike Fig. only the absolute value of the separation 
is considered here. A Kolmogorov-Smirnov test gives a 
probability p « 10“^ of both samples being drawn from 
the same distribution. As above, the finite resolution of 
our spectra censors the distribution of peaks with small 
separations, but the typical separations observed exceed 
our resolution limit by a large factor, suggesting that we 
are able to resolve the intrinsic distribution of peak sep¬ 
arations for each sample. Furthermore, our use of the 
same instrument setup and line-analysis algorithms for 
both the LAE and LBG samples allows us to demonstrate 
an instrinsic difference between the two populations. 

In contrast to the results of Yamada et al.| (2012a), 
we find no significant correlation between the Lya peak 
width and peak separation among multi-peaked pro¬ 
files for our LAE sample, and only a moderate asso¬ 
ciation among the LBGs (Pearson correlation r « 0.3, 
p « 0.05). However, these parameters are highly corre¬ 
lated (although with significant scatter) in the combined 
LBG-I-LAE sample, for which a Pearson correlation test 
yields p « 10“"^ (r « 0.3). 

Here we have defined the peak separation by the ve¬ 
locity difference between the peaks of the two most dom¬ 
inant components of the Lya emission line, rather than 
the flux-weighted centroids of those components. Be¬ 
cause the typical Lya emission lines have redshifted pri¬ 
mary peaks with asymmetric extended red tails (partic¬ 
ularly for the LBGs; Fig. [^, using flux centroids or sym¬ 
metric Gaussian fits to dSine the peak separation would 
lead to 1) a larger median peak separation for both sam¬ 
ples, and 2) a larger difference between the LAE and 
LBG distributions. 

For comparison, the distribution of nebular lin e wi dths 
for both the LA Es a nd LBGs is given in Fig. [T^ As 
discussed in Sec. |3.2[ the nebular line widths trace the 
velocity dispersion of the stars and H II regions within 
the galaxy. The larger nebular line widths of LBGs with 
respect to LAEs are therefore likely to reflect the larger 
masses of continuum-bright galaxies with respect to faint 
LAEs. Further discussion of the mass-dependence of the 
interstellar and circumgalactic properties of these galax¬ 
ies is reserved for Sec. [5l 

In contrast, the velocity distribution of emergent Lya 
flux depends most sensitively on the velocity distribu¬ 
tion and the optical depth of the scattering H I gas. 
The observed spectral morphologies therefore suggest 
that LBGs are associated with higher velocity and/or 



'^neb (fc™' « ^ ) 


Fig. 14. — Distribution of nebular line widths for LAEs (grey) 
and KBSS LBGs (orange) for the subset of objects with MOSFIRE 
rest-optical spectra. Dashed lines denote the average line width for 
each sample. The instrumental resolution of 35 km s~^ has been 
removed from each value; objects consistent with this resolution 
have been assigned = 0. A typical LBG has a nebular line 

width ~2x that of a typical LAE. 

more optically-thick outflows than their LAE analogs. 
It is thus evident that the observational characteris¬ 
tics that govern a galaxy’s selection by the narrowband- 
excess and/or continuum-color techniques (namely, the 
observed Lya and continuum emission) must be linked 
to properties of the outflowing gas - whether directly 
(by the coupling of Lya photons to interstellar and cir¬ 
cumgalactic gas) or indirectly (through the star-forming 
regions that produce Lya, far-UV, and ionizing photons 
and presumably provide the energy/momentum for large- 
scale galactic outflows). In the remainder of this paper, 
we separate these effects of gas velocity and optical depth 
on the emitted spectrum by considering the absorption 
and emission profiles in composite galaxy spectra. 


4. SIGNATURES OF OUTFLOWS IN STACKED SPECTRA 
4.1. Absorption signatures of outflows 

While the resonance of the Lya transition makes it a 
sensitive tracer of gas in and around high-redshift galax¬ 
ies, this strong coupling can obscure the details of the 
physical processes driving Lya absorption and emission. 
In particular, we noted in Sec. that Lya transmission 
is sensitive to a variety of factors, including the optical 
depth, covering fraction, dust content, and kinematics of 
the gas distribution. Because they produce degenerate 
effects on the Lya profile, Lya emission alone is insuffi¬ 
cient to fully characterize the physical conditions among 
LAEs or high-redshift galaxies generally. 

The measurement of absorption lines in the rest-UV 
continuum spectra of galaxie s is thus a crucial to ol for 
disentangling these effects. Kunth et al. (1998) ana¬ 


lyzed the spectra of nearby Lya-emittmg galaxies, find¬ 
ing blueshifted metal absorption features from enriched, 
outflowing gas. This study noted the correspondence of 
the Lya emission and metal absorption in tracing the 
velocity structure and porosity of the interstellar gas. 
At higher redshifts, the analysis of the continuum spec¬ 
tra of typical star-forming galaxies began with spectra 
of the gravitationally-lensed z ~ 2.7 galaxy MS 1512- 
cB58 (Pettini et al.l 2000 2002) and was extended to 
non-lensed galaxies thron g stacking the spectra of many 
(^1000) z ~ 3 LBGs by Shapley et al.| (|2003|). These 
studies utilized the high S/iN of the lensed or stacked 
spectra to extract metal abundances, ion-specific cover¬ 
ing fractions, outflow velocities, and systemic redshifts 
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for the galaxies while placing these properties in the 
context of their sta r-formation rates, m asses, and Lya 
emission properties. Steidel et al. (2010) presented both 
interstellar absorption line features and Lya emission in 
the context of a unihed kinematic model for gas outflows 
in LBGs (Sec. 5 of that paper). For faint LAEs such as 
those in our sample, constraints on the metal content and 
outflow velocities are especially interesting because they 
reveal the extent to which past and ongoing star forma¬ 
tion are already having an effect on the chemistry and 
kinematics of these particularly young, low-mass galax¬ 
ies. 

LAEs are by selection generally faint in the contin¬ 
uum, so abs orption measurements of comparable fidelity 
to those of Shapley et al. (2003) present a significant 
challege, parti cularly tor ob.iects as fa int as those con¬ 
sidered here. Hashimoto et al. (2013) studied the ab¬ 
sorption profiles in a stack ot tour lAEs with systemic 
redshift measurements to measure outflow velocities, and 


Shibuya et al. (2014) conducted a similar analysis for four 
individual bright lAEs, finding typical outflow velocities 
fabs ~ 100 — 200 km s“^. As noted in Sec.Jl] how¬ 
ever, the objects in these studies have B-banclmagni- 
tudes (me) ~ 24 (L ~ L*), meaning that they easily fall 
within typical selection criteria for continuum-selected 
star forming galaxies. While useful for studying the cor¬ 
relation between outflow velocity Uabs and WLya, these 
samples are clearly far removed from the faint LAEs con¬ 
sidered in this paper, with median continuum magnitude 
(ttib) = 26.8 (L ^ O.IL*, ~10x fainter than these previ¬ 
ous samples). 

Eor such faint objects, large samples are needed to ob¬ 
tain even a low-fidelity measurement of metal absorp¬ 
tion. To make such a measurement, we stack the rest-UV 
spectra of all 318 spectroscopic KBSS-Lya LAEs using a 
rest frame defined by Zgys.Lya = gLva.pe ak ~ 200 km s“^ 
based on the redshift calibration of Sec. |3.2[ As many of 
our objects have multiple LRIS spectra of their rest-UV 
continuum, we include all spectra for each object such 
that each spectrum receives equal weightin g, e ffectively 


weighting each object by total exposure timej^ Eurther- 
more, two spectra with individually-significant continua 
(possibly AGN) were removed from the sample in order 
to ensure that they did not dominate the stack; all the 
remaining spectra have S/N <§; 1 per pixel in the contin¬ 
uum individually. Lastly, the portion of each spectrum 
corresponding to A « 5577 ± 15A in the observed frame 
was omitted from the stack in order to remove contami¬ 
nation from a bright sky line. In total, 422 LAE spectra 
were stacked, producing a final spectrum with an effec¬ 
tive exposure time of ~675 hours (all from the Keck 1 
10m telescope). 

An error spectrum was then generated for the stack 
via a bootstrapping procedure. 1000 iterations were per¬ 
formed in which 422 spectra were drawn (with replace¬ 
ment) from our sample and stacked to make a set of 1000 
randomized spectral stacks. The standard deviation of 
values at each pixel was then used to define the error 
spectrum of the true stacked dataj^ Erom this error 


The results of the continuum stacking are insensitive to 
whether the objects are weighted equally, weighted by number of 
spectra, or weighted by total exposure time. 

The error spectrum thus includes the variance in continuum 


TABLE 4 

LAE Continuum Absorption Features 


Ion 

Alab" 

(A) 

f b 
Jose 

ITion" 

(A) 


(10* s-1) 

Absorption 

features 





Si IP 

1190.416 

0.575 

0.53 ± 

0.15 

— 

Si IP 

1193.290 

0.277 

0.24 ± 

0.11 

- 

Si II 

1260.422 

1.22 

0.53 ± 

0.10 

— 

O I 

1302.169 

0.0520 

0.43 ± 

O.lQf 

- 

Si II 

1304.370 

0.0928 

0.17 ± 

0.08f 

— 

C II 

1334.532 

0.129 

0.32 ± 

0.11 

- 

Si IV 

1393.760 

0.513 

0.30 ± 

0.12 

— 

Si IV 

1402.773 

0.255 

0.21 ± 

0.08 

- 

Emission features 





Si II* 

1264.738 

— 

— 


30.4 

Si II* 

1265.002 

— 

— 


4.73 

O I* 

1304.858 

— 

— 


2.03 

O I* 

1306.029 

- 

— 


0.676 

Si II* 

1309.276 

— 

— 


6.23 

C II* 

1335.663 

— 

— 


0.476 

C II* 

1335.708 

- 

- 


2.88 


^ Vacuum wavelength of transition 
Oscillator strength from the NIST Atomic Spectra Database 
(www.nist.gov/pml/data/asd.cfm) 

Eq uiva lent width of absorption in stacked LAE spectrum 

(Fig.[I^ 

^ Einstein A-coefficients from the NIST Atomic Spectra 
Database 

® The Si II AA 1190,1193 doublet is in the continuum blueward 
of Lya, which suffers from lower S/N and a high density of 
absorption and emission lines. As such, these features are only 
used here to constrain the saturation of Si II A1260. 

^ The O I A1302 and Si II A1304 absorption lines are partially 
blended, so the division between them is somewhat uncertain. 


spectrum, we estimate that the stack has a median S/N 
^ 10 per pixel in the range 1230A < Arest 1420A, 
the range of interest for the absorption lines discussed 
below 

Measurements were made of the absorption of si x m etal 
line^ corresponding to four spectral regions (Eig.f^ Ta¬ 
ble Q. Eor each spectral region, the continuum lew! was 
independently estimated using the local median value, 
but the continuum estimates among all four regions are 
consistent to < 2% in A/a units. The measured contin¬ 
uum was also checked for consistency with the broadband 
photometry. At the median wavelength of the absorption 
lines (Aobs ^ 5000A), the stacked continuum flux density 
is 7.5 X 10“® Jy (tosooo = 26.72), comparable to the 
mean Lya-subtracted broadband B/G flux (7.9 x 10“® 
Jy; (tobb) = 26.66) measured photometrically among 
the 318 objects in the spectral stack. Absorption mea¬ 
surements were then made for the low-ionization lines 
Si II A1260, O I A1302, Si II A1304, and C II A1334, as 
well as the higher-ionization doublet Si IV AA1393,1402. 


level at each pixel, so it is a conservative estimate of the error in 
the mean spectrum. 

Unfortunately, we are unable to obtain robust measurements 
of the C IV doublet 1549,155lA, generally one of the strongest 
absorption lines in high-redshift galaxies and a key indicator of 
AGN activity when seen in emission. Because we favor the res¬ 
olution of the 600-line grism over the increased spectral range of 
the 300 or 400-line grisms, the C IV doublet falls off the spectro¬ 
graph for many of our objects. Furthermore, the particular redshift 
range of our LAEs {z ~ 2.6) causes significant contamination near 
Arest ~ 1550A due to the bright sky line at AqI^q ~ 5577A. 
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Fig. 15. — Interstellar absorption signatures of metal lines in the stacked spectrum of all 318 spectroscopic LAEs. Vertical red lines 
are the rest-wavelengths of metal ions seen in absorption, while dashed blue lines denote the wavelengths of fine-structure and nebular 
emission features. A detailed list of these lines is given in Table Yellow shaded regions are the areas used to calculate the equivalent 
width and absorption-weighted velocity of each species; the absorption estimated for Si II A1304 is shaded green to differentiate it from the 
partially-blended O I A1302 absorption. Uncertainties for the equivalent width measurements are estimated via a bootstrapping procedure 
described in the text, and the bootstrapped error spectrum is displayed in grey at the bottom of each panel. 



Fig. 16.— Curve-of-growth analysis for the 4 detected Si II lines 
in the stacked LAE spectrum. The scaled equivalent width of ab¬ 
sorption for each transition i is given by VUi/A^/osc i'> which is pro¬ 
portional to the column density of the ion (and thus equal for all 
transitions) if the transition is unsaturated. If transitions are sat¬ 
urated, VFi/A?/Qsc I will decrease with increasing /og^ Errorbars 
are the l-a uncertainties on Wi from the bootstrap analysis. The 
dashed horizontal line is the best-fit constant value of VU/A^/osc 
for the three transitions with the lowest values of /osc (the least 
likely to be saturated). 


The equivalent width of absorption is estimated in each 
case by numerically integrating the set of contiguous pix¬ 
els with /norm < 1 nearest the rest-frame metal-line wave¬ 
length in the normalized stack. The O I A1302 and Si II 
A1304 absorption lines are partially blended, so the divi¬ 
sion between them was determined by visual inspection 
of the stacked spectrum. This division was cross-checked 
in bootstrapped composite spectra resampled at a range 
of wavelength scales. Uncertainties in the equivalent 
widths were computed by performing the same measure¬ 
ment for each line in each of the 1000 bootstrap spectra; 
the quoted error refers to the standard deviation of these 
bootstrap measurements. An absorption-weighted veloc¬ 
ity was then measured for each line via numerical Integra- 
tion of the line profile: v^hs = /At,[l “ fnorm{v)\vdv/Av. 
While the equivalent widths in absorption of the two lines 
in the Si IV doublet were measured separately, they were 
constrained to have the same total absorption-weighted 
velocity. 

All the measured absorption lines have mean veloci¬ 
ties — 100 < {vahs^ ^ —200 km s“^ with absorption ex¬ 
tending to Uabs ~ —500 km s“^ (where negative veloci¬ 


ties indicate blueshifted absorpti on), in good agreement 
wit h the brighter LAE sam ples offHashimoto et al. (20131 
and IShibuya et al.| (|2014|) . We therefore hnd clear evi- 
dence for blueshifted metal absorption signatures in these 
L ~ O.IL* LAEs, a strong indicator of the presence of 
metal-enriched outflowing gas. 

Because the systemic redshifts are not precisely known 
for the LAEs in our stack, there may be some additional 
uncertainty in our measured velocities compared to those 
measured in samples with complete catalogs of systemic 
redshifts. However, our stacked spectra also show sev¬ 
eral hne-structure a nd n ebular emission features (dashed 
vertical lines in Fig. 151 that are more likely to trace the 
systemic galaxy redshifts; the correspondence of these 
peaks to the rest-frame wavelengths of these features fur¬ 
ther suggests that our corrected Lya line measurements 
are a good proxy for the redshifts in these LAEs. No¬ 
tably, we have not corrected the line strengths or veloc¬ 
ities for emission from nearby fine-structure transitions; 
the effects of this contamination likely causes an under¬ 
estimate of intrinsic absorption and overestimate of the 
velocity of Si II A1304 and C II A1334. 

Galaxy outflows may also be characterized by the cov¬ 
ering fraction of the UV continuum by absorbing ma¬ 
terials and its optical depth. While these properties 
are typically degenerate for individual transitions (ex¬ 
cept at much higher S/N and resolution), they may be 
constrained independently when multiple transitions of 
the same species are present. In particular, the expected 
ratio of the equivalent widths of Si II A1304 and Si II 
A1260 is IEi 26 o/kEi 304 = 12-3 when both transitions are 
optically thin. The measured ratio of these lines in the 
LAE composite spectrum is IEi 26 o/kbi 304 = 3.6, consis¬ 
tent with a saturated Si II AI260 transition, with the 
caveat that the measured Si II A1304 absorption has 
small S/N and is partially blended with the O I A1302 
line. To clarify this measurement, we consider the ab¬ 
sorption due to the Si II AAl 190,1193 doublet in the 
rest-EUV (Arest 1200A) spectrum. The high density of 
absorption and emission lines and lower continuum level 
of the spectrum shortward of Lya (including foreground 
Lya absorption) result in larger uncertainties in these 
line strengths, but both lines of the doublet are detected 
with high significance. 
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TABLE 5 

Lyo Spectroscopic Subsamples and Composite Lva Emission Profiles 


Subsample 

(IFLy„> 

{n) 

-^obj 

-^spec^ 

(7j.ed (km s !)*> 

fTblue (km s 1)'" 

Av (km s 

Lya redshifts^ 

IVLyo, > 60A 

iloA 

27.6 

158 

211 

115±5 

191±10 

447±7 

20 < Wi^yc, < 60A 

37A 

26.8 

160 

211 

116±5 

276±15 

468±7 

IVLya < 20Ad 

iiA 

26.6 

82 

94 

119±7 

353±37 

475±10 

KBSS LBGs 

27k 

24.4 

65 

65 

178±10 

226±8 

646±20 

Systemic redshifts 

MOSFIRE LAEs 

Aik 

26.9 

32 

57 

131±9 

234±14 

490±23 

KBSS LBGs 

27k 

24.4 

65 

65 

215±20 

246±11 

568±25 


^ Because some objects were observed multiple times, and all available spectra were combined in each continuum 
stack, Nspec > ^obj • 

^ The parameters crj.gj, trijiue, and Av are the velo c ity widths of the red peak, blue peak, and peak separation 
in the fit to the composite Lyo spectrum (Figs. [T7] fc |18| |. The peak widths reflect the subtraction in quadrature 
of the LRIS instrumental resolution of = 100 km s“^. 

^ Composite spectra constructed using Ly« redshifts, shifted by —200 km s“^ 

Objects observed on the basis of a narrowband excess and spectroscopically-detected Lya emission line, but 
which have VPLya smaller than the adopted threshold of WLya 20A 


The measured absorption equivalent widths (W) for 
each of the four Si II transitions are given in Table 
along with the corresponding vacuum wavelengths (A) 
and oscillator strengths (/osc)- In the optically-thin 
regime, these quantities are related to the column den¬ 
sity of absorbing material by A^ion oc Wi/Xf fosc,i for 
each transition of index i. Because the physical col¬ 
umn density of a single ionic species must be single¬ 
valued, the quantity W/X'^f osc will likewise be constant 
for all ground-state transitions of the same ion unless one 
or more transitions have saturated. In this latter case, 
W /}? /osc will decrease with increasing /osc ■ 

The quantity osc is compared for the four Si II 

transitions in Fig.j^ The Si II A1190, A1193, and A1304 
transitions are consistent with a single value of W/X^fosc, 
but the Si II A1260 absorption falls significantly below 
the weighted average of the other three lines, indicat¬ 
ing saturation. While the other three transitions appear 
to be optically thin, both Si II A1193 and Si II A1304 
have nearby excited fine-structure emission transitions 
that may bias our measurements toward lower equivalent 
widths, as discussed above. If the Si II A1304 line were 
saturated, however, the bottom of its absorption trough 
would reach a similar depth as the saturated Si II A1260 
transition with respect to the loca l con tinuum (as is seen 
in some LBG samples; see Sec. 4.3), which is clearly 
not the case in the LAE composite spectrum (Fig. 15). 
This suggests that any filling-in of the absorption by ex¬ 
cited fine-structure emission has a small effect on the 
total measured absorption. The analysis of the com¬ 
posite spectrum thus strongly suggests that ti 26 o > 1, 
ti 304 ^ Ij and Tugo ~ Tugg < 1 based on the consis¬ 
tency of the curve-of-growth analysis in Fig. 16 (where 
T\ is the central optical depth of the Si II transition at 
wavelength A). 

Given the saturation of Si II A1260, the maximum 
depth of the line profile with respect to the continuum 
indicates a covering fraction of Si Il-enriched outflowing 
gas fc ~ 0.3 among the LAEs. While multiple unam¬ 
biguous absorption features are not available for other 
low-ionization species in the composite, the consistency 


of the absorption profile depths and widths for the O I 
A1302 and C II A1334 with the saturated Si II A1260 
line suggests that those species are co-spatial with Si II 
and have a similar covering fraction. We therefore take 
the covering fraction of low-ionization gas outflows to be 
fc ~ 0.3 for the LAE sample, noting the possible ex¬ 
istence of narrow absorption components unresolved by 
our R ~ 1300 spectra. 

The two lines of the Si IV doublet have a measured 
ratio bFi 3 g 3 /VFi 402 = 1.4, suggesting that the doublet 
is also saturated (VFi 3 g 3 /IFi 402 = 2 on the linear part 
of the curve of growth). However, the weaker line has 
S/N ~ 2, and ^25% of the bootstrap spectra have 
l^i 393 /bFi 402 ^ 2, so the transition may in fact be opti¬ 
cally thin. The higher ionization potential of Si IV (33.5 
eV) requires radiation from hot stars and/or collisional 
ionization in T > 10‘^K gas, and therefore traces more 
strongly ionized gas in or around the galaxy. The depth 
of the two absorption lines thus suggests the presence 
of an ionized medium with an average covering fraction 
fc ~ 20% (or greater if the lines are optically thin). It is 
unclear whether this gas is co-spatial with the medium 
giving rise to the lower-ionization lines described previ¬ 
ously; this relationship is probed further by comparing 
the velocity profiles of their absorption below. 

It is important to note that all of the above constraints 
on the covering fraction of metal-enriched gas rely on 
the assumption that the metal-line photons scattered out 
of the galaxy spectrum are emitted away from our line 
of sight to the galaxy. In the limit that the scatter¬ 
ing medium is entirely co-spatial with the observed light 
distribution from the galaxy (and that the scattering is 
isotropic and photon-conserving), we would expect that 
the net emission and absorption from scattering would 
be zero. In reality, many photons will be scattered at 
large radii with respect to the light profile of the galaxy, 
particularly for transitions where excited fine-structure 
states provide an alternative escape pathway with low 
optical depth. Given the spatial compactness of the 
LAEs, however, this “emission-filling” may reduce the 
observed absorption, particularly for those transitions 
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without nearby excited fine-structure states such as Si IV 
AA1393,1402. This effect may partially explain the lower 
measured value of fc for the Si IV transitions, in addition 
to their unknown optical depth. 

The velocity shift of Lya emission shows a strong in¬ 
verse trend with Lya equivalent width in the collective 


popu lations of high-redshift LAEs and LBGs (Erb et al. 


20141, but there are many mechanisms that could un- 


deriy this relationship. If a change in outflow velocity is 
the driver, then we should expect the velocity profile of 
the metal absorption lines tracing the outflows to vary 
with W^ya as well. Unfortunately, the S/N in stacked 
subsamples is too low to probe this effect in individual 
lines. However, we can boost the S/N by combining the 
line profiles for multiple lines that we expect to trace the 
same gas. In particular, we can stack the low-ionization 
lines (corresponding to more neutral H I gas) and the 
high-ionization lines (corresponding to more highly ion¬ 
ized H II gas) and thereby extract enough signal to com¬ 
pare the variation of the gas velocity distribution with 
fULya. 

Spectra combined in this way are displayed in Fig. [17] 
for three subsamples divided according to Wnya, with 
median WLya ~ uA, 37A, and IlOA respectively. Be¬ 
cause these stacks include features at multiple wave¬ 
lengths, we subsampled them by a factor of two with 
respect to the original wavelength scale before stacking 
in order to maximize the velocity resolution of the stacks. 
Details of these subsamples are in Table 

Several key features are apparent in the bottom pan¬ 
els of Fig. [TtI Firstly, there is substantial absorption by 
low-ionization species in all three panels at the systemic 
redshift of the galaxy (v = 0), but there is little or no 
corresponding high-ionization absorption. At the nega¬ 
tive edge of the absorption profiles (v ^ —500 km s“^), 
however, the high-ionization and low-ionization absorp¬ 
tion are well aligned in each panel. This pattern may 
indicate that both ionized and neutral gas are present in 
the material most likely associated with an outflow, but 
that the interstellar medium is dominated by neutral H I 
and low-ionization metals at u = 0. As noted above, 
however, the Si IV transition is particularly sensitive to 
emission-filling, which will tend to systematically reduce 
the observed absorption at z; ~ 0. 

Secondly, if the absorption lines are assumed to be 
optically-thick, there is a slight variation in covering 
fraction as a function of WLyo,: fc ~ 30% for the 
WLya llA sample, fc ^ 20% for the intermediate 
sample, and fc ^ 25% for the highest-WLya sample. 
The decline in fc from WLya llA to Wuya ^ 37A 
is c onsistent with expe ctations for Lya-emitting galax- 
(2010) demonstrate that observed Lya- 
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les: 

emission is closely correlated with the covering fraction 
of the gas through which these photons escape. The in¬ 
crease in apparent fc between the (Wuya) = 37A and 
(Wuya) = IlOA samples may not be significant because 
the continuum emission is extremely faint in the latter 
stack; not only is the statistical error large, but small 
systematics in the spectroscopic background subtraction 
could be significant in the stack of these faint objects, 
possibly biasing the measured continuum level. How¬ 
ever, the consistency between the photometric and spec¬ 
troscopic continuum estimates suggests this effect should 


be quite small. In any case, it seems that the apparent 
depth of the interstellar absorption features correlates 
only weakly with Wuya for the faintest, highest-Wnya 
bin of galaxies. 

Thirdly, there is some dependence in the maximum 
absorption velocity as a function of Wnya- While the 
absorption extends to z;abs ~ —500 km s“^ in all three 
panels, the Wnya ~ llA stack shows absorption in Si IV 
extending to higher velocities at low optical depths. Con¬ 
versely, the most-blueshifted edge of absorption in the 


Wuya IlOA stack extends no further than Uabs —400 
km s“^ H It is important to note that the systematic ve¬ 
locities used here are derived from a constant correction 
to the Lya redshift, and thus the evolution of the Lya 
offset with Wbya could theoretically mimic an evolution 
in absorption velocity with Wbya- However, despite pre¬ 
vious measurements of the significantly smaller Lya off¬ 
sets among LAEs with respect to LBGs, we find that the 
peak of the Lya profile has a consistent 200 km s“^ red¬ 
shift with respect to systemic at our observed resolution 
for both LBGs and LAEs (Sec. |3.2[). Furthermore, no 


difference with W^ya is observed m the most-redshifted 
edge of the low-ionization absorption (presumably cor¬ 
responding to neutral interstellar gas at the systemic 
redshift) across the three panels of Fig. such a shift 
would be expected from a Wbya-dependent estimator of 
the systemic redshift. Lastly, the total velocity width of 
absorption decreases significantly with increasing Wbya 
across the three panels, suggesting the presence of a real 
anti-correlation of Lya equivalent width with outflow ve¬ 
locity among these faint LAEs. 

Given that Wuya is likely to be anti-correlated with 
star-format ion rate and stellar mass among these LAEs 
(as discussed in Sec. 3.1), we interpret the above trends 
as a positive correlation between star-formation rate and 
outflow velocity. A similar analysis of t he K BSS LBG 
sample (Fig. 18 discussed further in Sec. 4.3) suggests a 
continuation of this trend to higher star-formation rates. 


4.2. Lya emission signatures of outflows 

The composite Lya e miss ion profiles for each subsam¬ 
ple (top panels of Fig. 17) show evidence for changes 
in outflow velocity with fVLya as well. Each Lya pro¬ 
file exhibits a narrow peak redward of the systemic red¬ 
shift (by construction, as they are stacked on the basis 
of their peak Lya redshifts) as well as an extended blue 
wing, which broadens with decreasing WLya (Fig. 17). 
To quantify this effect, a two-component Gaussian fit was 
performed to each stacked Lya profile with the following 
form: 


/A(^^)=^bluee-("+^"/2)/2"blue + 

, ( 6 ) 

such that the fit profile consists of two Gaussian peaks of 
arbitrary height and width, but with equal and opposite 
shifts with respect to systemic (±Az;/2). The model in¬ 
tentionally evokes an idealized outflow scenario in which 
the average LAE is surrounded by outflowing gas in both 

Note that the maximum absorption velocity is insensitive to 
any bias in the continuum level due to systematics in the back¬ 
ground subtraction. 
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Fig. 17.— Lya emission and interstellar absorption signatures in stacked subsets of spectra grouped by their photometric Lyo; equivalent 
width, Details of each subsample are given in Table Top panels: Lya emission spectra (black), and profile fit (green) for each 

subsample. The solid green line is the full fit, the dotted green line is the fit continuum level, and the dashed green line is the blueshifted 
Gaussian component of the fit. The fits do not extend to n > +500 km s“^, as this portion of the spectrum was omitted during fitting. 
Dashed and dot-dashed vertical lines denote n = 0 and v = —500 km s“^, respectively. Note that the specific flux in pjy is measured in 
the observed frame. Bottom panels: Red curves are the average absorption profiles of the low-ionization species Si II A1260 and C II A1334, 
while blue curves are the absorption profiles of the more highly ionized Si IV (A1393 and A1402). O I A1302 and Si II A1304 are omitted 
due to blending. As above, dashed and dot-dashed vertical lines denote = 0 and v = —500 km s~^, while the dashed horizontal line 
denotes the local continuum level. 




Fig. 18.— The composite Lyo: emission and low/high-ionization absorption line profiles (as in Fig. |17[ l for the KBSS LBG sample. 
Left: The composite spectrum is constructed using the nebular (systemic) redshift measurements. Ri ght: Ltie composite is constructed by 
stacking spectra according to their Lya emission peak, in the same way as the LAE samples of Fig. |l7| As in Fig. the composite Lya 
profile is distorted with respect to the composite constructed using systemic redshift measurements. However, the composite absorption 
profile is not significantly affected by the choice of nebular redshifts or calibrated Lyo redshifts. 
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directions along the line of sight. In all three stacks, the 
red side of the Lya profile diverges strongly from a Gaus¬ 
sian shape at large velocities, forming a broad red tail 
similar to that seen at u ^ 0. In order to maintain the 
simplicity of the model, the spectrum at u > -1-500 km 
s“^ was omitted from the fitting. 

The resulting fits are displayed in Fig. with fit pa¬ 
rameters given in Table Both the width of the blue 
component and the separation of the blue and red peaks 
increases with decreasing WLya- While the Lya peak 
width and separation may be increased by raising the 
H I column density, such a scenario would not predict 
a corresponding increase in the velocity width of metal 
absorption transitions, which have significantly smaller 
optical depths than Lya. The similarity of the extended, 
blueshifted Lya emission to the metal absorption pro¬ 
file therefore indicates that emission from outflowing gas 
likely populates the wings of the Lya emission profile. 
That is, we suggest that the extended Lya emission 
observed at large redshifts and blueshifts with respect 
to systemic depends sensitively on the outflow velocity, 
rather than merely the optical depth of the scattering 
medium. Models for this behavior, including the depen¬ 
dence of Lya transmissivity on both the amplitude and 
vel ocity gradient of outflowin g gas, are discussed in detail 
bylSteidel et al.|(|2010| |2011|). 

Radiative transfer simulations of Lya transmission 
suggest that these processes can generate double-peaked 
emission profiles with asymmetric offs ets with respect to 
the sy stemic redshift. In particular, |Verhamme et 
(2006) and others find that a spherical shell of gas ex¬ 
panding with velocity t4xp can in many cases produce a 
blue peak at v = — 14xp and a red peak at v = 2I4xp- 
The stacked Lya profiles in Fig. (or Fig.[I^ discussed 
below) are insensitive to the peak of the blue component 
and cannot discriminate between the two models (that 
is, whether the offset of the blue peak is equal to that of 
the red peak or half that of the red peak). However, the 
same qualitative relations hold for either model: the sep¬ 
aration of the two peaks and the breadth of the extended 
blue and red components increases with increasing con¬ 
tinuum luminosity. 

4.3. Comparison to continuum-bright galaxies 

While the trend of the inferred outflow properties with 
Wnya is weak among the LAEs in this sample, we can 
probe a much larger range in galaxy properties by ex¬ 
tending the sample to include continuum-bright galax¬ 
ies. In particular, the composite Lya profile of the 
KBSS LBGs (Fig. 18) displays a continuation of several 
trends identified among the LAE subsamples in Fig 17 
The KBSS LBGs all have systemic redshift measure¬ 
ments, so we can compare the results of stacking them 
by either their nebular redshifts (Fig. left panel) or 
their Lya-derived redshifts (right paneT^as in Fig. fl7|). 
While the detailed profile depends on which redshittis 
employed, the LBG composites have broader and more 
widely-separated peaks than any of the LAE subsam¬ 
ples, and the velocity extent of the broad tails of the 
Lya emission again corresponds closely to the velocity 
extent of blueshifted absorption from the higher- and 
lower-ionization species. N ote that the absorption pro¬ 
files (bottom panels of Fig. 18) are insensitive to which 
redshift indicator is used: tfe^ Si II A1260, C II A1334, 



Fig. 19.— Composite UV continuum spectra for stacks of KBSS- 
Lya LAEs and two comparison samples of LBGs. Details of the 
LBG samples are given in Table|^ The LAE stack (black, bottom) 
is the same as that in Fig. |15| and the bootstrapped error spec¬ 
trum is again shown in grey (primarily visible near the wavelength 
of the Lya line). The KBSS LBG stack (orange, middle) has red¬ 
shifts measured via nebular emission lines. The “LyC LBG” stack 
(purple, top) has redshifts measured by a calibrated combination 
of Lya emission and metal absorption redshifts. All three compos¬ 
ite spectra have R ~ 1300 for 1250A< A < 1400A. The challenge 
of measuring outflow signatures in faint LAEs is clear: not only is 
the continuum emission ~10x fainter than typical LBGs, bu t th e 
underlying absorption signatures are significantly weaker (Fig.|20|l. 


and Si IV doublet absorption centroids each change by 
less than ~25 km s“^ when stacked according to the neb¬ 
ular or calibrated Lya redshifts. This consistency sug¬ 
gests that our measurements of interstellar absorption 
velocities in the LAE composite are unaffected by the 
scarcity of systemic redshifts among the LAE sample. 

In Figs. 19-21 we compare the stacked spectrum of 
all 318 LAEs in this sample with two stacks of LBGs: 
the KBSS LBG (Steidel et al.|2014) sample, and another 
sample from a survey tor Lyrnan-continuum (LyC) emis¬ 
sion in z ^ 3 galaxies (C. Steidel et al., in prep.). A 
comparison of these two samples to the KBSS-Lya LAE 
sample is given in Tabl e |6l _ 

As discussed in Secs. |2.2| & |3.2[ the KBSS LBG sam¬ 
ple consists of 65 galaxies that each exhibit Lya emis¬ 
sion in their UV spectra and have accurate systemic red¬ 
shifts measured from rest-frame optical emission lines. 
These rest-UV spectra were observed using Keck/LRIS- 
B with the 600-line grism (the same used for the KBSS- 
Lya LAEs), providing a spectral resolution R ~ 1300, 
or (Tinst ~ 100 km s“^. There is significantly stronger 
metal absorpt ion in each of the spectral lines shown 
in Figs. 15 & [T^ and the absorption profiles are both 
broader (^tending to u « —1000 km s“^) and deeper 
(<0.5x the continuum value) than the corresponding 
features in the LAE stack. The excited fine-structure 
and nebular emission lines (Fig. 20), however, are well- 
matched to the corresponding LAe features, further val¬ 
idating the redshift offset employed for the KBSS-Lya 
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Fig. 20.— Composite rest-FUV continuum spectra (as in Fig. |l9[ | normalized to a common /1325 continuum level (1 250 A< A < 1400A). 
As in the preceding figures, the bootstrapped LAE error spectrum is in grey. The absorption lines displayed in Fig. |15| are here marked 
by vertical dashed lines. The two LEG samples (orange, purple) are highly consistent with each other and exhibit clear contrasts with the 
LAE sample, including significantly broader and deeper absorption profiles. 



Fig. 21 .— Composite rest-EUV continuum spectra (as in Fig. |l9| l normalized to a common /1125 continuum level (1050A< A < 1200A). 
Note that the normalizations differ from Fig. |20[ the continuum ratio / 1325//1125 ^S'Ch stack is shown here by the horizontal lines 
labeled “FUV” on the far left of the figure. The KBSS LEG composite is omitted due to the low S/N in this region, but the continuum 
ratio for that spectrum is marked with the orange horizontal line. As in the preceding figures, the bootstrapped LAE error spectrum is in 
grey. Insterstellar absorption lines are marked with vertical dashed lines, color-coded by atomic species and labeled from above. Horizontal 
dashed lines connect transitions of the same ion. Vertical solid lines labeled from below mark features with a photospheric component. 
Significant blending is visible due to both the absorption features and the many excited-state emission features throughout this region of 
the spectrum (not marked). 


TABLE 6 

Comparison of LAE and LEG Continuum and Absorption Properties 


Sample 

(^> 

mi450‘‘ 

(mag) 

Mi450’^ 

(mag) 

'^Lya'^ 

(^) 

-^obj 

(A) 

(km s“^) 

(A) 

■wcir 

(km s“^) 

(A) 

■“SilV® 
(km s“^) 


KBSS-Lya 

2.70 

27.0 

-18.3 

42.7 

318 

0.53 

-98 

0.45 

-185 

0.51 

-190 

0.3 

KBSS LBGs 

2.36 

24.4 

-20.7 

27.2 

65 

1.67 

-224 

1.35 

-201 

1.98 

-276 

0.6 

LyC LBGs 

3.07 

24.6 

-21.0 

21.5 

74 

1.50 

-172 

1.44 

-176 

2.41 

-288 

0.6 


^ 7711450 observed AB magnitude corresponding to the rest-frame A ~ 1450A UV continuum. At 2 : ~ 2.7 (KBSS-Lya and 

KBSS LBGs), this correponds to the G band, whereas the 7^ magnitude is used for the LyC LBGs at 2 ~ 3. For consistency with 
the other samples, the LAE 7711450 value is not Lya-subtracted, though the Lya emission dominates the broadband flux in many 
cases. 

^ Mi 45 o is the absolute magnitude corresponding to 7711450 the typical redshift of the sample {z). 

^ Because the LBG samples have no narrowband Lya images, all three values of WLy^, here were measured spectroscopically for 
consistency. Note that the spectroscopic LAE sample has a mean photometric equivalent width Wrya “ 85.6A. 

^ ^SilL ^CIL WsiTv are th e eq uivalent widths in absorption of Si II A1260, C II A1334, and Si IV A1393-I-A1402 (combined), 
respectively. See Tableland Fig. |15| 

^ Tfsiii, TJQii, and ugiiY are the absorption-weighted velocity of Si II A1260, C II A1334, and Si IV A1393-I-A1402. 

^ fc is the covering fraction of low-ionization gas implied by the depth of the saturated Si II A1260 absorption trough, which is 
consistent with the covering fraction implied by the depths of the O I A1304 and G II A1334 transitions in each composite spectrum. 
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stack. 

The LyC LBG sample consists of 74 galaxies, most of 
which do not have measured nebular redshifts; the red- 
shifts for this sample were estimated via a combination of 
their Lya emission and interstellar absorption redshifts 
with a cali bration based on th e results of ISteidel et al.l 
(2010) and Rakic et al. (2011). Furthermore, the com¬ 
posite spectrum was checked after stacking using stel¬ 
lar photospheric features. About 15% of these galaxies 
have redshifts derived solely from the Lya line and as¬ 
sume a 300 km s“^ offset (with respect to the Lya cen¬ 
troid) based on stacks of the subset of objects with neb¬ 
ular emission redshifts. Like the LAE and KBSS LBG 
comparison samples, only objects with spectroscopically- 
detected Lya emission were included in the sample. The 
correspondence of the absorption and emission signatures 
of the stack to those of the KBSS LBG stack again 
demonstrates the efficacy of the calibration based on 
UV features. Being at higher redshifts than the other 
samples, these spectra were observed using a combina¬ 
tion of the LRIS-R 600-line grating (covering the spec¬ 
trum Arest ^ 1200A) and the LRIS-B 400-line grism (for 
Arest ^ 1200A). This lower resolution was employed to 
maximize the S/N of the spectra shortward of Lya (and 
the Lyman limit), where the observed spectrum becomes 
particularly faint. These spectra also have signihcantly 
larger integration times compared to the KBSS LBGs 
and KBSS-Lya LAEs (8—10 hours vs. 1.5—2 hours). 
Additionally, the higher redshifts of this sample (z ^ 3) 
with repect to the KBSS LBG sample (z ~ 2.3) cause 
the lOOOA < Arest < 1200A continuum to lie at observed 
wavelengths with significantly higher LRIS-B through¬ 
put, further boosting the S/N of the LyC LBG compos¬ 
ite spectrum at these wavelengths with respect to that of 
the KBSS LBGs. This stack is thus particularly useful 
for comparing the absorption due to Ly/3 and other EUV 
transitions among the LBG and LAE samples (Fig. 21). 

The absorption profiles of the two LBG composite 
spectra are quite similar (Fig. 20). The velocity ex¬ 
tent of the absorption is almost identical in both sam¬ 
ples (umax ~ —1000 km s“^), though the LyC LBG stack 
shows slightly deeper absorption across all the lines. This 
difference may signify a difference in the galaxy proper¬ 
ties: at z ^ 3, the LyC LBGs have a brighter aver¬ 
age rest-UV luminosity than those at z ^ 2.3 (which 
have similar apparent magnitudes). The total equiva¬ 
lent width in interstellar absorption (Table [ 6 ), which is 
insensitive to the spectral resolution, is sinmar for the 
LyC and KBSS LBG stacks. 

While the blending of Si II A1304 with O I A1302 in- 
hibits a direct meas urement of the Si II optical depth, 
Shapley et al. (20031) demonstrated that the Si II A1260 
transition is typically saturated in similar samples of 
bright LBGs, and the absorption of the Si II A1190, 
A1I93, and A1526 transitions in the LBG stacks are con¬ 
sistent with this interpretation. Furthermore, the con¬ 
sistency of the absorption profile depths among the low- 
ionization species in the LBG stacks suggests a gas cover¬ 
ing fraction fc ~ 0.6 for enriched, neutral/low-ionization 
material for the two LBG samples. The observed line ra¬ 
tio of the Si I V doublet is IFi 333 |W ^ Am ~ 2 for both LBG 


stacks (as in Shapley et al. 


thin absorption and /c ^ 0.4 tor the ionized H II gas. 


2003[), suggesting optically- 


The difference in absorption line strength between the 
LBG and LAE samples can likewise be seen in the Ly/3 
transition (Fig.[2T|). While the LBG composite is almost 
completely absorbed at the Ly/3 line, the LAE spectrum 
has significantly higher transmission (> 50%). The con¬ 
tinuum ratio / 1325//1125 (displayed as horizontal lines on 
the left side of Fig. 21) also varies significantly between 
the LBG and LAE samples: the LBG composites have 
a flux ratio / 1325//1125 ~ 1-6, while the LAE compos¬ 
ite has / 1325//1125 ~ 1-3. The flux ratio for the z ^ 3 
LyC LBG composite is uncertain because the spectra at 
Arest ^ I 2 OOA and Arest ^ 1200A fall on different detec¬ 
tors (LRIS-R and LRIS-B) within the spectrograph, so it 
is not entirely clear how much of the observed variation 
in continuum absorption is due to physical differences in 
the LAE and LBG galaxy populations, rather than the 
changing tran smissivity of the IGM with redshift (e.g.. 


Madau 1995). Both effects are likely to be present, as 
the U V continuum level blueward of Lya is particularly 
sensitive to the dust content and metallicity of the ISM 
(the latter due to metal-line blanketing) and the age of 
the stellar population, all of which are observed to be 
lower among LAEs than LBGs. 

Through comparison of the LAE and LBG stacks, it 
seems clear that the qualitative links between L ya emis¬ 
sion and and interstellar absorption discussed in|Shapl^ 


et al. (2003) persist down to far fainter continuum lu¬ 
minosities tnan those probed by samples of LBGs. The 
outflow velocity (represented by either the absorption- 
weighted mean or the highest-velocity edge of absorp¬ 
tion) decreases with increasing W\,ya across the LBG and 
LAE subsamples, and the covering fraction of gas drops 
steadily with increasing to at least the intermediate 
subsample of LAEs (WLya 37A). The offset of the Lya 
emission centroid likewise decreases from the LBG sam¬ 
ples (-1-300 km s“^) to the LAE stacks (-1-200 km s“^). 
While further observations and modeling are required to 
determine the physics of star-formation driven outflows, 
we suggest that our results provide an important link be¬ 
tween the physical properties of these galaxies and those 
of the gas in and around them. We provide some initial 
interpretation of these results in the context of previous 
work below. 


5. DISCUSSION 

As described above, the UV continuum luminosities of 
these galaxies are closely linked to their star-formation 
rates, and our brighter samples of galaxies (e.g., LBGs 
and low-WLya LAEs) are associated with higher veloc¬ 
ity dispersions and thus larger dynamical masses. As¬ 
suming that the star-format ion rate (SFR) and stellar 
mass (M*) are roughly proportional to the UV luminos¬ 
ity 0-6., the luminosity corresponding to M 1450 in Ta¬ 
bled, the masses and star formation rates of the KBSS 
LBGs likely exceed those of the LAEs by a typical fac¬ 
tor SFRlbg/SFRlae A7*_lbg/A7*,lae 10 . Simi- 
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differ by a factor ^1.5—2 between the LAE and LEG 
samples: the Lya emission cent roid ( i^la e ~ 200 km 
s“^, vlbg ~ 300 km s“^; Figs. 17 & 18), the absorp¬ 
tion centroid of low-ionization lines (wlae ~ —100 km 
s“^, 'Clbg ~ —200 km s“^; Table and the maximum 
velocity of blueshifted absorption (ulae ~ 500 km s“^, 
800 km s“^; Figs. 17 & 18). This factor of 


fLBG 


A.5—2x in velocity over a lOx range in mass and star- 
formation rate is broadly consistent with the scaling seen 
in samples of galaxies at at low and intermediate redsh ifts 
( MartinpOOS Rupke et al.||2005 Weiner et al.ll 


which i^outflow is measured to scale as This 

positive scaling stands in contrast to the lack of outflow 
velocity variation with eithe r star-format ion ra te or ve¬ 
locity dispersion observed by Steidel et al. (2010). In that 
paper, the authors note that their LEG samples have 
st ar-formation r ates similar to the highe st-SFR galaxies 
of Martin (2005) and Rupke et al. (2005), where the rela- 
tionship between SFR and ^outflow flattens off. Our LAEs 
are likely to have similar star-format ion rates to the low- 
SFR galaxies of the low-redshift samples (SFR < 10 Mq 
yr“^), so our results may indicate that the relationship 
between star-format ion rate and outflow velocity (i.e., 
positive scaling at low SFR, saturation at SFR > IOMq 
yr“^) is similar at low and high redshifts. 

Analogously, the doubling of Uneb between the LAE 
and LEG samples here suggests that outflow velocity 
grows in a roughly linear (or marginally sub-linear) fash¬ 
ion with the stellar dynamics of the galaxies. Similar 
scaling is observed in the lower-redshift samples noted 
above. Furthermore, if t he escape speed sc ales as Vesc ~ 
5 —buneb (as discussed in Weiner et al. 2009), then the es¬ 
cape velocities of the LAEs'7iri~i4^7'«2DlL This velocity 
is of order the typical absorption-weighted centroid veloc¬ 
ity of the outflows (and significantly less than the maxi¬ 
mum velocity of absorption), suggesting that a significant 
fraction of outflowing gas is permanantly ejected from 


Weiner et al. 

2009 

). More detailed measurements of the 

masses and s 

;ar-tormation rates ot individual LAEs are 


required to probe the scaling of these feedback processes 
within the sample of low-mass, high-redshift galaxies. 

It may seem paradoxical that the centroid or width of 
the Lya line is linked to outflow velocity and/or the star- 
formation rate in a galaxy while the Lya peak velocity at 
fixed resolution shows no dependence on continuum lumi¬ 


nosity. In particular. Figs 
pling between the position 


seem to show a decou- 
narrow Lya peak and 


the breadth of the extended Lya emission tails to pos¬ 
itive and negative velocities. We suggest that previous 
studies of Lya emission actually show a similar decou¬ 
pling between the Lya peak and the width of the overall 
velocity distribution of Lya flux . The radi ative-transfer 


model s employed by Kulas et al. (2012) and Ghonis et al. 
(2013) seem to show that the dominant red peak of Lya 


transmission through an expanding spherical shell of gas 
does not sensitively depend on gas velocity (at least for 
Ugas ^ 300 km s“^), but the distribution of flux in sec¬ 
ondary and tertiary red and blue peaks becomes more ex¬ 
tended with increasing gas velocity. These narrow emis¬ 
sion line features are observed in higher-resolution spec¬ 
tra of bright a nd n earby galaxies ( [Verham me et al.|2008 


ful for modeling the effects of radiative transfer and gas 
velocity directly, but they are not feasibly observable in 
L ^ O.IL* galaxies at z > 2. We therefore suggest that 
the observed uniformity of the Lya peak offset is due to 
the absorbing presence of some neutral gas at u = 0 in all 
of these star-forming galaxies (which sets the minimum 
velocity shift for an escaping Lya photon), whereas the 
emission at larger (positive or negative) velocities reflects 
those photons scattered off of the high-velocity outflow¬ 
ing gas (which depends more strongly on galaxy prop¬ 
erties). The broad wings of the Lya emission line (e.g., 
Fig. are therefore likely to be populated by scatter¬ 
ing oiFof outflowing gas, and thus are sensitive to the gas 
velocity and show a velocity distribution similar to that 
of the extended metal absorption. 

However, the relatively low measured covering frac¬ 
tions of these outflows and the consistency of the Lya 
peak offset may suggest that the majority of Lya scat¬ 
terings in LAEs occur without incurring any single large 
kinematic “kick” from an interaction with the highest- 
velocity gas. This tendency would be likely in a sce¬ 
nario where galaxy-scale outflows are highly anisotropic 
and thus create an orientation-dependent model for Lya 


outflows is 

highly favored by simulations (e.g., iHopkins 

et al.| 

2012) and lower-redshift observations (e.g., Kornei 

et al. 

wn 

at z ^ 1), both of which suggest that out- 


disks, but it is not clear whether the anisot ropy of the 
outflows and that of Lya emissivity coincide. |Verhamme| 


et al. (2012) suggest that the Lya flux and equivalent 
width are maximized when a galaxy is observed face- 
on. However this model would naively predict a positive 
correlation between Lya velocity offset and equivalent 
width (rather than the anti- or absence of correlation 
observed) because both the Lya emission and the line- 
of-sight outflow velocity would be maximized for face-on 
galaxy disks. 

Local star-forming galaxies can help probe the geome¬ 
try of Lya escape. In particular, the Lya Refere nce Sam¬ 
ple (LARS; Hayes et al.|2013 Ostlin et al.|2014 ) consists 
of 14 local galaxies with Hb'r imaging and spectroscopy 
of their Lya emission. This sample was selected to pro¬ 
vide a set of local analogs to high-redshift LAEs that can 
be ob served at much higher spatia l and spectral resolu¬ 
tion. iRivera-Thorsen et al. (2015) use HST/COS spec¬ 


tra to analyze Lya emission along with absorption due 
to Si II and Si IV along the line of sight to the brightest 
star-forming regions within these galaxies, finding sev- 
eral results broadly co nsistent with those of this paper. 
IRivera-Thorsen et al.| use four spectrally-resolved Si II 
transitions to show that the typical absorbing medium is 
optically thick, with a covering fraction that depends on 
both position and velocity. As in this paper, they find 
that higher local (and global) Lya transmissivity is asso¬ 
ciated with low covering fractions of Si II, the Si Il/Si IV 
absorption ratio is high at u = 0, and the resolved Lya 
emission peak shows little dependence on outflow veloc¬ 
ity. 13 The authors interpret these results in light of a 


Chonis et al.||2013 Martin et al.||2(J15| ) and may be use- 


In contrast, |Henry et aL| (|2015ll argue that the covering 
fraction of neutral gas does not dorninate Lya escape in low- 
redshift, lower-mass galaxies, basing their claim on the appearance 
of Lyman-series absorption lines with zero residual transmission 
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model where Lya escapes through an expanding medium 
that is both porous and non-uniform in velocity, with in¬ 
homogeneities potentially driven by Rayleigh-Taylo r in¬ 
stabilities in expanding bubbles. Hayes et al. (2014) find 
that the Lya emission is spatially extended with respect 
to the FUV continuum emission (in agreement with high- 
redshift studies); the local Lya/Ha ratio is higher at 
larger galacto-centric distances. In some edge-on disks, 
Lya/Ha is seen to peak above and below the disk plane, 
but in most LARS galaxies the ratio appears to vary with 
no simple depe ndence on galaxy geometry or structure. 


Furthermore, Law et al. 


(2009) find that most galaxies 


at 2 ~ 2 — 3 are not dynamically dominated by rota¬ 
tion and thus are less likely than low-redshift galaxies 
to generate coherent, collimated outflows. The young, 
faint populations of LAEs considered in this paper may 
be even less likely to maintain such structural coherence. 
Therefore, the KBSS-Lya LAEs may be best modeled 
by irregular morphologies in which Lya photons escape 
through a patchy ISM with minimal planar symmetry. 
Such a model, in which Lya photons escape through 
“holes” in the neutral gas distribution, would similarly 
explain the pos sible association betw een Lya and LyC 
emissivity (e.g., Mostardi et al.|[2M3 ). 

LBGs and other bright star-forming galaxies can have 
extremely low Lya and LyC escape fractions even when 
low-resolution spectroscopy suggests that the gas cover¬ 
ing fractions are fc < I; in many cases, there may be ex¬ 
tremely narrow, unresolved absorption components cov¬ 
ering the line of sight that are still capable of absorbing 
all of the incid ent LyC flux (e.g.,th e g alaxy Q0000-D6, 
as dis cussed by |Shapley et al .||2003| and |Ciallongo et al.| 


20021. Opening up lines of sight through the neutral ISM 
that are transparent to both Lya and LyC photons may 
be easier for low-mass galaxies. 

Other mechanisms of Lya escape could lead to decou¬ 
pling between the Lya and LyC emissivity of galaxies. 
While the scattering of Lya photons is driven by H I, the 
eventual absorption and destruction of these photons is 
dominated by dust. If the LAEs contain little interstellar 
dust and have velocity fields capable of quickly scatter¬ 
ing Lya photons out of resonance, then a large fraction 
of Lya photons may escape before encountering a dust 
particle (without requiring a “hole” in the H I distribu¬ 
tion) . If Lya escape is dominated by this channel, then 
the Lya and LyC escape fractions may be less strongly 
correlated, given that LyC photons are absorbed by H I 
rather than merely scattered. The apparent breakdown 
of the correlations of galaxy properties with Wnya for the 
highest-WLya bin may indicate the transition to such a 
regime where dust formation and mixing into the ISM 
has not fully progressed. The faintest LAEs at 2 ^ 2 — 3 
may be most analogous to those observed at z ~ 6, where 
the extremely blue UV slopes suggest that escaping pho¬ 
tons encounter very little dust while exiting the galaxy. 
More robust fitting of the LAE spectral energy distribu¬ 
tions, including deep NIR images from MOSFIRE, could 
strongly constrain the dust content of these young galax¬ 
ies. Similarly, deep NIR spectra may constrain the metal- 


in cos spectra of Lyg-emitting Green Pea galaxies (GPs; |Car-| 
|damone et al.||2009f. Ho weve r, the KBSS LAEs show signihcant 
transmission at Lyp (Fig. HR, suggesting that the gaseous envi¬ 
ronments of GPs are qualitatively different from those of 2 : ~ 2 — 3 
LAEs despite their similarly low masses and compact sizes. 


licity (and thus the chemistry) of their ISM. 

Another possibility is that the decoupling between Lya 
emission and galaxy properties in this sample at the high¬ 
est values of WLya may result from the contribution of 
the QSO ionizing radiation field to their Lya flux. Eor 
fluorescently-illuminated systems, the total Lya flux is 
proportional to the surface area of optically-thick H I 
and the strength of the local QSO-generated i onizing field 
( Cantalupo et al.|2005 Kollnieier et al.|20I0 ), with mini- 
mal dependence on the kinematic or chemical properties 
of the galaxies. Our forthcoming analysis of the HST 
NIR and optical images of these fields will shed more 
light on the physical sizes, morphologies, and stellar com¬ 
ponents of the Lya-emitting regions of these galaxies. 

6. GONCLUSIONS 

This paper has presented the rest-frame UV spectral 
properties of Lya-emitters near hyperluminous QSOs at 
z ~ 2.7 from the KBSS-Lya survey. These LAEs are 
significantly fainter (L ^ O.IL*) than previous spec¬ 
troscopic samples of high-redshift star-forming galaxies, 
which enables us to extend observations of interstellar 
enrichment and kinematics to a new galaxy regime. Our 
primary results are as follows: 

• The faint LAEs have Lya escape fractions /esc,Lya ~ 
30% (based on the observed Lya to Ha flux ratio), sig¬ 
nificantly higher than previous observations of b righter 
star-forming galaxies at high-redshift (Sec. JSL This 
result is consistent with the conclusions of ISteidel et al.l 

dMIl). 


• Like continuum-bright galaxies at similar redshifts, 
LAEs produce Lya emission that is significantly shifted 
from the systemic redshift. While the flux-weighted 
Lya emission centroid varies with galaxy luminosity 
across our combined sample of LAEs and LBCs, the 
Lya peak (when observed with R « 1300) displays 
a highly consistent shift of ^^200 km s“^ irrespective 
of continuum luminosity. This offset may be used to 
estimate redshifts of diverse Lya-emitting galaxies ob¬ 
served at comparable resolution whe n re st-frame opti¬ 
cal spectroscopy is unavailable (Sec. 3.2). 


• Faint LAEs display a range of diverse Lya morpholo¬ 
gies and are associated with narrower Lya peaks with 
smaller peak separations than those of LBCs. These 
results are consistent with smaller velocities and/or 
column densities of neutral gas in the circumgalac- 
tic medium of LAEs compared to LBCs. While most 
of the Lya emission lines have a characteristic red- 
dominant morphology, there is a minority ( 13%, 4 ob - 
jects) with blueshifted primary peaks (Secs. 3.2-3.41. 


• By combining the UV continuum spectra of all 318 
LAEs in our spectroscopic sample, we find absorp¬ 
tion signatures of several metal ions that have been 
shown to trace outflowing material in more luminous 
star-forming galaxies. We believe these are the first 
measurements of metal-enriched outflowing gas in such 
faint high-redshift galaxies. We measure absorption- 
weighted velocities of ~ 100 — 200 km s“^ and max¬ 
imum velocities of ~ 500 km s“^, both significantly 
smaller than those seen in comparison samples of LBGs 
(Sec.|4]). 
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• Using multiple transitions of Si II, we infer the satura¬ 
tion of the Si II A1260 transition, thereby constraining 
the covering fraction of low-ionization material in front 
of the LAEs to fc ~ 0.3. This covering fraction is sig¬ 
nificantly smaller than that observed in LBGs, with 
fc ^ 0.5 — 0.7 in the same ions. This low covering frac¬ 
tion of obscuring gas likely drives the observed high es¬ 
cape fractions of L yg a nd Ly-continuum photons from 
faint galaxies (Sec. 4.1). 


• We find that outflow velocity - whether probed by 
maximum velocity of metal absorption, absorption- 
weighted mean velocity, or Lya emission morphology 
- increases with continuum luminosity (decreases with 
Lya equivalent width) among LAEs and with respect 
to LBGs. This result suggests that these gas outflows 
are tied to physical properties of their associated galax¬ 
ies - e.g., the mass and star-formation rate, for which 
the nebular velocity dispersion and continuum lumi¬ 
nosity are proxies. In particular, our results are con¬ 
sistent with weak dependence of outflow velocity on 
star-formation rate (^outflow SFR°-^®), and we find 
that the outflow velocities (up to 500 km s“^) are large 
compared to the probable ~200 km s“^ escape veloci¬ 
ties of the LAEs (Sec. [^. 


In sum, these observations demonstrate that stellar 
feedback is an essential ingredient in realistic simulations 
of even the youngest, least massive galaxies at high red- 
shift. In particular, these feedback processes are key to 
understanding the escape of line and continuum photons 
(particularly ionizing photons) from faint galaxies and 
their effect on their surrounding IGM. Further work is 


necessary to analyze the apparent absence of correlations 
between Wi^ya and covering fraction of ISM absorption 
at the highest equivalent widths. This effect may be 
the result of a qualitative transition in the properties of 
galaxies to those with minimal enrichment and dust, and 
could also point to the effect of the QSO in picking out 
galaxies independently of their star-formation rates. 

In each of these topics, further progress will be driven 
by continued near-IR observations of faint galaxy pop¬ 
ulations. Rest-frame optical spectra (in combination 
with observations in the rest-UV) are the only way to 
effectively constrain the star-formation rates, metal en¬ 
richment, and gas kinematics of these newborn galaxies, 
thereby establishing the relationships between their stel¬ 
lar and gaseous properties and characterizing evolution 
of the baryon cycle over cosmic time. 
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